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A Geometrical and Mechanical Method to determine 
Lorentz Transformations 


SUMÁRIO 


Apresenta-se um método geométrico de obtenção 
das Transformações de Lorentz numa representação 
gráfica em que os eixos coordenados e homólogos dos 
dois referenciais se fazem coincidir, Utilizam-se 
coordenadas cartesianas ortogonais. 

O método é susceptível de uma fácil mecanização 
tendo em vista o traçado automático de trajectórias 
nos diferentes planos coordenados. 

Sob outra forma e numa diferente aproximação, 
o método foi apresentado pela primeira vez numa 
comunicação feita em 6 de Maio à classe de Ciências, 
da Academia das Ciências de Lisboa, sob a designa- 
ção «Método geométrico de obtenção das Transfor- 
mações de Lorentz. Propriedades do método e suas 
aplicações». 

Num próximo número da «Técnica» será publi- 
cado um artigo, em português, no qual se procederá 
a uma aplicação do método à resolução de determi- 
nados problemas básicos da electrodinâmica dos cor- 
pos em movimento. 


|. PRESENTATION OF THE METHOD 


Two inertial frames of reference are considered: 
s(x, y, Z, t) and S(X, Y, Z, T), both being rectilinear 
orthogonal Cartesian coordinates with homologous 
axes parallel to each other. The velocity of s, relative 
to S, is v in the direction of the positive X-axis, The 
origins of s and S coincide at the time t=T=O. 

In this case the Lorentz Transformation is: 
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M. ABREU FARO 
Electrodynamics Centre (|. N, |. C.) 
Technical University of Lisbon (1.S.T.) 


SUMMARY 


A geometrical method to obtain Lorentz Trans- 
formations in a graphical representation where the 
axes of the two frames of reference are made to 
coincide, is presented. Orthogonal coordinates are 
used and the method can be mechanized with a view 
to an automatic drawing of trajectories on the diffe- 
rent coordinate planes. 

Under a different form and with a different 
approach, the method was firts presented to the 
Science Class of Academia das Ciências de Lisbon, 
on the 6th of May 1976, under the title «Geometri- 
cal Method to obtain Lorentz Transformations. 
Properties of the Method and its Applications». 

In one of its next issues Tecnica will publish 
an article in portuguese where the method will be 
shown as applied to the resolution of certain basic 
problems of the electrodynamics of moving media. 


X=y (x+ 8 ct) 


cT=y (ct+ px) (1) 
Y=y &=z 
where 
A v | 
S — —— and = reed 
C y | — ps 


and c is the velocity of light in vacuo, 
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In this method we shall make the t, T-axes coin- 
cide with y, Y-axes, see Fig. 1. 


Fig. 1 


Hlustrates the essencial features of the method Through 
directions + y/2 and points P, and Q, a bicontinuous 
mapping is established from [x, t] space into [X, 1] space 


In order to obtain symmetrical expressions with 
respect to space and time, we shall consider c = 1. 
After this, (ct) and (CT) are reduced to t and T 
although having the physical dimension of a length, 
as obvious. 

Having assumed c = 1, the transformations (1) 
take the form 


X=y(x+ pt) 
(1 a) 


T=y (t+ px) 
where t and T have the dimension of a length, as 
stated, but measuring the time in metric units of 
length. 

The length and time units are the same in both 
frames of reference, In all figures v/c = 0,6. 

Our purpose is to present a geometrical method 
of determining the transformations 


[x, t] = [X, T] (2a) 


in a graphical representation where the homologous 
axes of the two frames of reference are made to 
coincide. 


Once this has been done the transformations 


[x yVI= IX, 3] 
(2 Db) 


[x, 2] = [X, Z) 


become immediate considering that by (1) Y=y and 
Z = 2. 


The method in question can be easily mechanized 
and thus contribute, we think, towards a simple and 
quick handling of the Lorentz Transformations spe- 
cially suited for examining the correlations between 
space-time coordinates (x, y, 2, t) and (X, Y, Z, T) 
and associated phenomenology. 

Beginning with Fig. 1 we shall first present the 
method to determine the transformations (2 a); after 
which the theory will be given. 

In order to obtain [X, T] from [x, t] we proceed 
in the following way: 

We consider the directions r,, R, r, and R, 
forming angles 


A (lr, t)=—y/2 
«+ (Ra, T) = + y/2 


«A (Tx, x) =—y/2 
A (R,X)=+W/2 


in accordance with the positive senses adopted and 


v 
where sin y=> — = ph. 
e 


Now, from the point [x, t] we draw a straight 
line with an r, direction; the point this straight line 
intersects the t-axis is P,; similarly Q, is the point 
an r, direction intersects the x-axis, 

Next we take P, and Q, as origins of R, and R, 
directions and these lines cross exactly at the [X, T] 
point. 

Conversely [x, t] could be obtained from [X, T] 
by using the same method. 


H. THEORY AND ASSOCIATED PROPERTIES 


ES gd fig 
The theory of the process is a very simple one. 
As it can be easily verified 


pis 
tan y/2 = a (3) 
WY 


From (la), the following Lorentz Transforma- 
tions can also be easily obtained: 


(= —+» q—1 ' 
[= a )=[ T+] (1) 


[es ft) fem] ue iate-s)] 


which means that: 
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i) Under the transformation (1a) the image of 
the straight line r;, ( — y/2) which intersects x = o 
at t=r7, is the straight line R, (+y/2) which inter- 
sects X=o at T=17,. 

ii) Similarly, the image of the straight line 
r, (—vy/2) which intersects t=o at x =1 is the 
straight line R, (+ y/2) which intersects T = o at 
X =. 

So essentially the method is this: 

If [x, t] belongs to the straight lines r, and r,, 
[X, T] belongs to their transforms which are R, 
and R, defined by P, and Q, which are common 
points of intersection at the t, T-axes and x, X-axes 
of the straight lines r, and R, and the straight lines 
r, and R,, respectively. 


Some interesting properties are associated with 
the established correspondence. 

Let us start by verifying that angle y/2 is the 
only one which satisfies the visualized symmetry of 
the transformed straight lines. 


£ 


at 
Let m and M be the slopes m = — 
AX 


aT 
M= AX of a straight line and its image respectively, 
A 


then from (1a) results; 


E 
1 + mê 
If we make 
m=z—-M 


an equation is obtained the roots of which are 


A | 
ai Ram ndo = tan y/2 
141 
M,=—m = = cotan y/2 = Ee (6) 


once (3) is considered. Obviously m, and M, corres- 
pond to r, and KR, and m, and M, to rr, and R.. 

It is also easy to demonstrate that m, and m. 
are the only directions for which a straight line and 
its image intersect respectively at the common points 
P, and Q,, on the t, T-axes and x, X-axes, 

These facts are a joint proof of the uniqueness 
of y/2. 


Another important property becomes apparent 
from Fig. 1. 

Let «a» be a point in s (x, t) and «<A» its image 
in S (X, T), then: 
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i) If the point «b» belongs to the straight line 
(a, r,], then «<B> is a point of its transform [A, R,). 

Considering the directions of r,, r, R, and R, 
an immediate conclusion is that 


ab = AB 


ii) In the same way, if «c» belongs to [a, r,] 
«O» belongs to [A, R,] and 


ac = AU 
— 
Hence, vector ad has in the axes f[a, (r,, r,)] the 


same components as vector AD in the axes TA, 
(Ro Rol. 

As a corollary, and since the transforms of the 
origins «o» and «O» are coincident, in our graphical 


—, 
representation the components of a vector oa in axes 


[o, (r; r.)] coincide with the components of OA im 
axes [O, (R, R,)] where «A» is the image of «a» 
and conversely. 

The first conclusion is particularly interesting 
to study the evolution of the points in the vicinity 
of «a» and «A». 


Finaly Fig. 2 makes it clear that major problems 


Fig. & 


The directions intervening are always rz and re and R, and 
Re. P and Q are the images of p and q. Once known the 
hybrid coordinates (x, T) (X, t), (X, x) or (T, t), [x, t] 
and [X, T] are easily obtained on the basis of this geometry 


can be solved taking as a basis y angles symmetri- 
cally arranged in relation to the coordinate axes, 
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i) Nameley, the following chains clearly stand 
out in Fig. 2 


Fig. 3 


Directions x = const., t= const. X = const. T = const. and 


their transforms, Velscity straight lines + c are transformed 
into vclocity straight lines + c 


xy>QakTIEAPaAlX tt] 
[A; €) 2 q =t [x €] 2 p + [XT] 


and solve the problems: 

Knowing [x, T] or [X, t], determine [X, T] and 
IF 

ii) Also, considering that the ordinates of P and 
pare (T+ t)/2 and the abscissae of Q and q are 
(X + x)/2 the following problems are also solved: 

Knowing [X, x] or [T, t], determine [X, T] 
and [x, tJ. 

These geometrical relations are general and easy 
to demonstrate, for simplification purposes [x, t] 
and [X, T] are both in the first quadrant, 

Through them it is possible to solve any problem 
using only the method presented, 

However, for the solution of certain problems, 
it is easier to draw the transformed straight lines 
from determined directions such as: 


x=const.”, t=const.”, X=const,.”, T=const.' 


* Reference [1]. 
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This was done in Fig. 3. The procedure will be taken 
up soon in an article where it will be duly used for 
a simple and quick resolution of certain concrete 
problems. 


Il. PHYSICAL INTERPRETATION OF THE METHOD 


in [1]* a different approach was used to obtain 
the method. Though it will not be repeated here, 
certain concepts then used will, nevertheless, be 
taken up. 

From now on the method which has appeared 
to be purely mathematical will also possess a physi- 
cal substratum in perfect agreement with the prin- 
ciples of Relativity. To this purpose let us consider 
the transformations (1) which, once differentiated, 
lead to 


AX = y ( 4x + pe dt) 
(8) 
CAT = y (cat + pAx) 
Let us have 
AX AX 
sa too dá 
(9) 
Ee — id 
a o í: sx 


As it is obvious U and u represent velocities in 
the frames S (X, Y, Z, T) and s (x, y, z, t), res- 
pectively. 

D and d represent that which from now on, will 
be designated «desynchronization»: variation of the 
coordinate time, T or t, by unit of length of the 
coordinate space, X or x. 

From (8) and (9) it is easy to obtain the well 
known law of composition of colinear velocities: 


uv 
a | + uv,c? Ea) 


and the equivalent relation 


peotrvo 


qr (10b) 


which at this stage does not vet possess any physical 
significance. 


From (8), (9) and (104) there result the follo- 
wing facts, always present in the deduction of the 
Lorentz Transformations. 

Let us fix our attention, along the time, on a 
point x, of the x axis. 

We shall have, successively: 
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x = const.* = x, 4x=0o, u=0, U=v 


which is interpreted in S (X, Y, Z, T) in the following 
way: 

s (x,y, Z, t) is animated with a velocity v rela- 
tively to S (X, Y, Z, T). 

If we fix X = X,, we would obtain u = —v 


Let us know make a different experiment. 

Let us concentrate our attention, along the x axis, 
on an instant t, of the coordinate time of s (x, y, Z, t). 

We shall have, successively, 


t=const.º'=t, At=o, d=o, D=v/€ 


which is interpreted in S (X, Y, Z, T) in the following 
way: 
s (x, y, Z, t) is animated with a desynchroniza- 
tion v/c2 relatively to S (X, Y, Z, T). 

If we fix T='T, we would obtain d =— v/c2, 

This symmetry is, like in the case of the velo- 
cities, in perfect agreement with the principles of 
Relativity. 


This being the case, to say that s (x, y, Z, t) 
is animated with a velocity v in relation to S (X, Y, 
Z, T) is as legitimate and equivalent a statement as 
to say that s (x, y, Z, t) is animated with desyn- 
chronization v/c? in relation to S (8, Y, Z, T). 

The first statement is present in the Galilei 
Transformation whereas the second is not: they are, 
nevertheless, equivalent. 


An immediat result of this is an assymetry in 
this transformation, that is: singular treatement for 
variable time which then assumed absolute character. 


For such ideas to materialize it will become ne- 
cessary to conceive experiments whereby the measu- 
rement of desynchronization will become accessible. 
In fact, a measurement that can be effected as natn- 
rally as the relative velocity v. 

For the measurement of velocity v let us take 
on x, a source of light which emits along the time 
a flash of a physically infinitesimal duration in rela- 
tion to the periods of repetition. 

In S (8, Y, Z, T) the successions 


“Tou "Tau Ke ... 


Koor Ho Mas =.» 


can be observed, where 
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No — Xo Xo  Xoy 


= DOS Ne - v (11) 
To — Too To — Tg 


In order to measure the desynchronization it will 
suffice to synchronize n sources along the x axis in 
such a way that they will be triggered simulta- 
neously at instant t,, producing n flashes at the 
time. 


InS(X,Y,Z, 'T) we will observe the successions 


O 


é TA ' To , Ros esenerensas 


where 


Tos— Too Too — Tor e 
Xo— Xowa No — Xo 


The question becomes even more interesting if 
we take, as previously, c = 1. 

In this case velocities and desynchronization are 
adimensional physical values and the relative velocity 
and relative desynchronization are measured by the 
same figure «v>». 


cce. V/C (12) 


However, attention must be drawn to the fact 
that whereas in the first experiment it would be 
possible for a material particle to travel along the X 
axis, accompanyng at each T,, instant the X,, image, 
in the second experiment this is not possible, To a 
desynchronization 


D=v/Cº <c 
there would correspond a velocity 
U=ce/v>c 


But recording clocks in X,, Xoa Xu... would 
register T,, To Ta... and would reveal a history 
with a chronology where each event could never be 
accompanied by the same particle, even if this par- 
ticle were a photon, 

The experiment, however, is feasible exactly as 
it was conceived since there is no difficulty in pro- 
voking in s (x, y, Z, t) simultaneous flashes of n 
sources nor is there any intrínsic difficulty in the 
conception of recording clocks. 

Ultimately this means that the straight lines 
X=>X(T) in the first experiment are time-like 
«whereas in the seccond they are space-like. 


Bearing all this in mind and assuming c=1, 
the expressions (10a) and (10b) are converted into 
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u+- vw 
— 13 
= l+-uv E 
d--ov 
am (13b 
o l+dv , 


With regard to the first of these expressions, 
we can formulate the following problem: 

Determine the velocity —u, ins (x,y, Z t) which 
can be transformed into a velocity u, in S (X, Y, 
o EM, 

It is obvious that the problem can be solved. 

From (13a) this velocity can be determined by 
imposing 
[uj=-u, 


[U1=+u, 


Hence 


-U + v 


| = (14) 


1 — uv 

Since the desynchronization composition law is 
the same, the following problem is immediately 
solved: 

Determine a desynchronization —d, in s (x, Yy, Z, 
t) which will lead to a desynchronization d, in S (X, 
das E) 

Furthermore, the equation which solves this pro- 
blem is the equation 


-di+v 


di= 
: I— dyv 


(15) 


formally identical with (14) and resulting from (13 b) 
where we imposed 


[D] = +d, [d]J=-—d, 


Since (14) and (15) are second degree equations 
it follows that the correct solution is in (14) 


u<v 
and in (15) 


d,<v 
As it would be expected, the other solutions 


U >v 
or 


d>v 


represent, by their inverse, the solution for the pro- 
blem of desynchronizations in terms of velocities 
and the resolution of the problem of the velocities 
in terms of desynchronizations. 


* 
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Solving equations (14) and (15) we obtain 


q—l 
E 


Since straight lines of the form 


P=MX+HT, 
and 
T=M(X—X,) 


were used, the slopes of the straight lines of velo- 
city, + u, are given by 


À agp A cet? = 2cotan y/á 
M, == By 


+M, = + 


whereas the slopes of the straight lines of desynchro- 
nization + d, are given by 


em =xtÊ se md 
By 


which coincides with the expressions (6) 


* 


In this way the straight lines r, and R, repre- 
sent the result of the following and easily conceivable 
experiment: 


Let us establish in s (x, y, Z, t) a chronology 
of events with —d, desynchronization (it will suffice 
to tell the history toward the past); S (X, Y, Z, T) 
will then show a chronology of events ordered inver- 
sely in accordance with d, desynchronization. 


The straight lines r, and R, represent, in their 
turn, the result of the following experiments: 


Let us put in s (x, y, Z, t) a particle travelling 
along the x axis at a —u, velocity; this particle will 
be observed in S (X, Y, Z, T) travelling at a u, velo- 
city along the X axis. 


These experiments are graphically represented in 
Fig. 4. 


This subject will soon be taken up in a publica- 
tion. 


For the time being we have just tried to give 
physical contents to each one of the pairs 


lr; Ri] and [ry RJ] 
which represent in themselves a feasible experiment 
in regard to the measurable concept of desynchro- 


nization and the more trivial and corresponding con- 
cept of velocity. 
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E e E E q De A E a 7 E + a E E 


Fig. 4 


To the sequences [a, 5, 6, 7] and [A, 5, 6, 7] there corres- 

pond, respectively, velocities — u; and + us. To the sequen- 

ces [a, 1, 2, 3] and TA, 1, 2, 3] there correspond, respectively 
desynchronizations — d; and + d; 


IV. MECHANIZATION OF THE METHOD 


Let us take two triangular plates, of aperture y, 
overlapping and with their vertices sliding along the 
coordinate axes: if one of the intersections (r,, TF.) 
is made to follow point [x, t], evolution [X, T] is 
immediately obtained, Conversely, if intersections 
[R, R,] are made to follow [X, 'T], evolution [X, t] 
is reproduced. 

In Fig. 5 the method is clarified in relation to 
two possible situations. 

From now on it becomes easy to conceive ana- 
logical and automatic processes to visualize relati- 
vistic effects and arrive at an automatic drawing of 

T=T(X) 


TY=Y(X) Z=Z (X) 


knowing 
=t (x) 


y=y (x) Z=z (x) 


V. FINAL CONCLUSION 


The method presented represents a geometrical 
method to obtain Lorentz Transformations. 
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MEXA NOR 
0a NS 

ANO 

SOS AN 


Fig. 5 


From overlapping plates of aperture y with their vertices 
sliding along the coordinate axes, results a mechanical pro- 
cedure establishing the relation between [x, t] and [X, T] 


Since the coordinate axes are made to coincide, 
the concept is different from that of the Minkowsky, 
Synge, Amar, Loedel and Brehme diagrams. 

In fact, these authors obtain the transformations 
on the basis of the determination, on non-coinciding 
axes, of the coordinates X and T or x and t. 

In our case, the axes coincide and the phenomena 
are presented without distortion, exactly as they can 
be «seen» by the observers of both of the inertial 
frames. 

We believe the facility with which the method 
can be mechanized will serve an easy handling of the 
Lorentz Transformations and, consequently, an appre- 
hension of the relations in the space-time continuum. 

The method has already proved to be rather use- 
ful and quick in the resolution of problems connected 
with the propagation of electromagnetic waves such 
as the study of aberration of light, Doppler effect, 
etc. 

The method in question also allows an easy inter- 
pretation of infinitesimal transformations, namely 
that of the significance and criterium of the quasi- 
-validity of the Galilei Transformation, 
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Algumas reflexões sobre representações matriciais de grupos 
de transformações 


RESUMO 


Dado um grupo de transformações e um conjunto 
de funções linearmente independentes discute-se o 
modo como, tomando como base as funções referidas, 
deve ser definido um conjunto de matrizes constituin- 
do um grupo isomorfo do grupo de transformações. 


1 — Admitamos que, dado um conjunto 1, 4, 
Va, de funções e uma transformação A de um grupo 
de transformações, a matriz associada a A quando 
as funções são tomadas como base se define por (1). 


3 
AU = > u aki (1) 
k 


em que a,, é o elemento genérico da matriz indicada, 
matriz que designaremos por IT (A). 

Desenvolvendo a expressão anterior para o caso 
das 3 funções v,, vu, 4, vem 


Avw = ay % + ar ta + ag Yg (2) 

Ava — Ajo Vy + Ago Vo + Ago Lg 

Avs = ag k% + agito + ag Ya 
De acordo com esta convenção de escrita, o con- 
junto « deve ser entendido como uma matriz linha 


e a multiplicação dessa matriz pela matriz associada 
a A deve ser escrita sob a forma 


A (vo tg = (by vo Ya) aj Ap AjgN = 


dao, &go Agg 


às Ago gg 
(3) 
= (ar W + ao Vo Ag Yg aj2% + Ago do + 


+ Aga bg 


2,3 %1 + Aos vo + Agg Vo) 
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(Investigador do IAC. Prof. aux. equip. no IST) 
Centro de Quimica Fisica Molecular 

Complexo Interdisciplinar 

Instituto Superior Técnico 


ABSTRACT 


Given a transformation group and a set of li- 
nearly independent functions, the method by which 
a set of matrices which forms a group isomorphous 
with the transformation group must be defined is 
discussed. 


Se em vez da transformação A tivermos as 
transformações B e € será analogamente 


3 
Hj sb 
AA > E a (4) 
k 
3 
Cy = U 
f by E” 5) 
k 
B (4 vo va) = (4 Vo 43) bj bjo Dis 
ba boa Dos 
bg Doo Ds 
(6) 
= (by +Dbata+Da ts Dry + Dota + 
+ bao ig Dist + Dogvo + Dgs Yg) 
C (tvotg) = (Vy$>%3) Cm Cj Ci 
Coy Com Co 
Cy Co Cy 
(7) 


= (CyW+Ccarto+ Ci Cry + Cato + 
+csotg CigW + Cog% + Cas 3) 


Suponhamos que as transformações A, Be € 
estão ligadas pela relação 


C = BA (8) 
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Isto significa, de acordo com a convenção habi- 
tual, que aplicando a uma função a transformação 
A e ao resultado a transformação B se obtém o 
mesmo resultado final que se obtém aplicando direc- 
tamente a transformação C., 

Apliquemos à expressão (3) a transformação B. 
Entrando com a matriz de transformação T(B) defi- 
nida em (6) tem-se, identificando o resultado com (7) 


BA (ly Yo vg) = (ag ty + avo agita Ago 4 
+ asa ba + Aga Ya Aya ty + Aos + Ag ta) X 


x D;, Dj Di es 
Da, Das Da 


Day Das Das 
(9) 


= (Cut + Cata taty Crab + Cor ba + Cos Ya 
Cia W + Cos Po + Cos ba) 

em que 

Cut + Cata Cuta = (AY + Ata Ag Ya) Da 
+ (aja % + Aga to + Aga ba) Day 
+ Agt + Anda + Aba) Da 

Cab Coba Cab = (AY + Anta + Ag Ya) Dia 

+ (ao Y + Alado Aga Pa) Das 

+ (Bb + Basta + Ag ba) Da 

Cat + Canta + Cuba = (A 0 + Aya Vo + Ag Ya) Dis 
+ (aa + Ao Vo + Aga bg) Dog (10) 
+ (Ag t + Aagbo + Aga ba) Da 


ou seja 
Cuty + Cada + Cata = (a Di + aja Da + Ag Dai) 4 
+ Ay Di + Ago Day + Aga Da) Va 
+ (ag Di + Ago Day + Asa Days) Va 
City + Colo + Cata — (Ay Dia + Aja Das + Ass Dao) 
+ (Ag Dig + Dos Dis + Aa Dag) Va (11) 
+ (Ag Dio + Ago Dos + Ag Da) Ya 
Cat + Cava t Cab = (a Di + Ara Day + Ag Dag) by 
+ (aja Dia + Boo Dog + Ago Das) da 
+ (ag Dis + Aga Das + Agg Ds) Ya 


De modo genérico tem-se 


3 
o == 
ij 2 x e (12) 
o que mostra que 
IMC)=T(A) T(B) (13) 
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ou seja que a matriz associada a € se obtém multi- 
plicando as matrizes associadas a A e B, pela ordem 
A —B, em contraste com a ordem B— A do produto 
das transformações indicado em (8). 

Como é óbvio pode chegar-se a este mesmo re- 
sultado usando as expressões de definição (1), (4) e 
(5). Na verdade por aplicação da transformação A 
às funções 4; obtém-se, de acordo com a definição (1) 
funções &, dadas por 


3 
, Gi » ' E" (14) 


Se aplicarmos agora B às funções à, obtém-se 
funções x, dadas por 


E eme > + b (15) 
i ] j Ji 


Mas de acordo com (14) 


 - » ns ai 


Substituindo (16) em (15) vem 


| sm b 
, ea E “ Er j 


1 
= W 
BA A 4 K. a" ) (17) 
k j 
ou seja 
L — ! 
i 2 " E” (18) 
com 
Cc = b 
ki 2 a ji (19) 


de acordo com (12) e (13) 
Consideremos no entanto de novo a expressão (1) 


4 DMA (20) 


Se aplicarmos a transformação B vem 


i 


BAY = B (24 a) 
k k ki 


= B 
> A rs (21) 
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Se atendermos a (4) pode escrever-se 


E b 
A ', A (22) 


Substituindo (22) em (21) obtém-se 


aaa sa > - ', E E” 


Este resultado parece mostrar, ao contrário do 
resultado a que se chegou anteriormente que o efeito 
global da aplicação da transformação A seguida da 
aplicação da transformação B é traduzido por uma 
matriz de transformação que se pode obter a partir 
do produto das matrizes de transformação associadas 
a Be A, por esta ordem, 


2 — Vejamos se é possível compreender a causa 
desta aparente contradição, Ao usar expressão (1), 
ela deve ser entendida de uma forma simbólica, e 
exprimindo o facto de que por aplicação da transfor- 
mação A à função “; do conjunto v,, “, vw, se obtém 
uma nova função, combinação linear das funções v,, 
Vos Ugo Os coeficientes dessas combinações lineares 
podem ser agrupados constituindo uma matriz, que 
será, por definição, a matriz associada à transforma- 
ção A que se designa por T(A). Tal como está apre- 
sentada, a definição (1) implica que as funções vi 
se devem entender como formando uma matriz linha. 
Por seu turno a matriz T(A) deverá ser uma ma- 
triz quadrada, de ordem 3 (igual ao número de fun- 
ções “; que constituem a base de representações). 
Para se poder efectuar a multiplicação das matrizes 
da forma habitual é necessário escrever primeiro a 
matriz das funções w e depois IT(A). Nessas condi- 
ções é lógico que se escreva, na verdade, o somatório 


do 2.º membro da expressão (1) sob a forma Sly a. 
pk ki 


Mas nesse caso, para manter a convenção, parece 
ser melhor, escrever, de modo simbólico 


fia ) ns (24) 


com a convenção que o primeiro membro significa 
a aplicação de A a y. Ter-se-ia então simbolicamente 


(Vi Va ba) À = (4 Va val TIA) (25) 


Se agora aplicarmos a transformação B deve 
usar-se a mesma convenção. Em termos matriciais 
já se viu (cfr. expressão 17) que aparece o produto 
das matrizes pela ordem A —B. Em termos de ope- 
radores, e na sequência da equação (25), deve escre- 
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ver-se a transformação B à direita, embora se deva 
admitir que a operação se efectua da «direita para 
a esquerda», isto é, aplicando primeiro A e depois B. 
D« qualquer modo escrever-se-á, simbolicamente 


(bi va va) AB = (VV Ya) T(A) TB) (26) 


Este resultado parece mostrar que do ponto de 
vista simbólico a aplicação da transformação A se- 
guida da aplicação da transformação B deve ser 
escrita sob a forma AB, colocada à direita das fun- 
ções 4, Essa transformação global (que pode con- 
tinuar a ser designada por C) teria uma matriz 
associada TI(C), produto das matrizes associadas a 
A e B, pela ordem AB, ou seja, a ordem inversa 
da aplicação das transformações, A expressão (26) 
escrever-se-ia então | 


WBC = (baby) PIC) (27) 
com 


que substituiria (8) embora continuando a significar 
o mesmo que (8), isto é, que €C é a transformação 
que aplicada a uma dada função tem o mesmo efeito 
que a aplicação de A seguida da aplicação de B. 


3 — Neste conjunto de raciocínios há contudo um 
ponto importante que importa esclarecer. Como se 
sabe, a matriz associada a uma transformação depen- 
de do conjunto de funções de base que se considera. 

Ora ao escrever a expressão (14) admitiu-se que 
ao aplicar B às funções 4, (que, recorde-se, tinham 
sido obtidas por aplicação de A às funções 4 ) a 
matriz de transfc.” «ção I(B) era a que tinha sido 
obtida quando us fuuções «w eram usadas como base 
(cfr. expressão (4) ). 

Atendendo às expressões (2) e ao significado 
das funções “ pode escrever-se 


é, = ay + Ayuda + Agby 
do = Ah + Aagta + Asgba (29) 
da = Ah + Aasvo + Bagda 


Quando se aplica B às funções (p, 6, 4,) tem-se 
Bo, = a, By + ayBy + ag Bb 


Bj, = a, By, + A, BY, + Aga Bl, (80) 
Bo, = a,By, + AnByo + AB 


| 


Mas de acordo com (4) 
By, a D,% + Dava + Dava 


Bys = Ds” je Dava + Dagta (31) 
By, o Dat + Dasta Ea Dagty 
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Substituindo (31) em (30) vem 


Bd, = ay (by + Dos Yo + Dar ba) + Bay (Dia by + Das Ga + 
+ Dag Ya) + Am (Dis Pi + Das Ya + Das Ya) 

B 6, = Alba by + Dada + Dar bs) + Aos (Dia + Das Ya + 
+ Dos ba) + ago (Dia by + Das bat Das Pa) (32) 

Bd, ag(Dya dy + Da va + Dar ba) + dg (Dia W + Dao Ga + 
+ Das va) + Aga (Dig Yi + Das da + Das Ya) 


ou de outro modo, 
3 3 4 
Bó = (A a box+ (2 a bo+ (2 a box 
3 3 3 
Bó — (2 aby+ (2 ab) + (2 ab x (83) 


3 3 3 
E W by 
» ' (A PM r 5) cr Ada (2 a RA 


Este não é ainda o resultado final que se pre- 
tende, uma vez que para determinar a matriz asso- 
ciada a B quando setomam como base as funções q, 
é necessário saber qual o resultado da aplicação de B 
às funções à, expresso sob a forma de combinações 
lineares das funções q&,. Para isso é evidentemente 
preciso exprimir as funções 4 em função das fun- 
ções q: 

Ora as equações (29) podem exprimir-se em lin- 
guagem matricial. Como é óbvio, deve escrever-se gy 
da € qn ASSim como v, v, e vw sob a forma de ma- 
trizes linha. 


(q, Pa Pg) = (4 W Ya) a,, a =” 
do do dos (34) 
a “yo 33 
Isto é 
(9, Po 03) = (My Ya Ya) TA) (35) 


Assim sendo, para exprimir os « em termos dos 
4 basta usar a matriz inversa de T(A) que represen- 
taremos por T-(A), Tem-se 


(vi Pa Ya) = (9, 04 P)L-MA) (36) 


Por outro lado, as expressões (33) podem escre- 
ver-se sob forma matricial 


B(q, Po Pa) = (4, Va Ya) 
ED A Dry Aya ED ds 


ED ds Ebay Ao Ear ds 


EDax do XD do Es By 
k k k 


(37) 


Isto é, por definição das matrizes T(B) e P(A) 
B(p, do Pa) = (W Va ty) POB) IA) (38) 
Substituindo agora (36) vem 
B(p, 4, 04) = (4, 4, JT A)T(B)LCA) (39) 
Este resultado é importante e mostra que quan- 
do se usam como funções de base as funções 4, à 
matriz que deve ser associada à transformação B não 
é T(B) (definida em (4) ) e que é a matriz associada 
à transformação B quando as funções |; são usadas 


como base, mas sim a matriz ["(B) dada, de acordo 
com (39), por 


P(B) =T+ (A) I(B) (A) (40) 
Nestas circunstâncias, em (9) não se deve usar 
a matriz T(B) mas sim a matriz [(B) de elemen- 
tos 
p=Fa, db 
a DER km 2mj (41) 
km 


sendo ix o elemento genérico da matriz T“(A) (que 
não é evidentemente, o inverso do elemento as). 
Daqui decorre que (12) deve ser substituída por 


8 
= "Dra. (42) 


ou seja, em vez de (13) deve escrever-se 


P(C) 


P(A)T(B) (43) 
Substituindo o valor de "(B) dado por (40) vem 
P(C) = IT(A) TI(A) T(B) I(A) 


ou seja 


P(C) 


r(B) T(A) (44) 
o que está de acordo com (8). 


Do mesmo modo, em (15), (17) e (19) os ele- 
mentos b,; devem ser substituídos por ESTE Nessas 
condições a expressão (19) deve ser substituída por 

Sa da Sa dá; (45) 


Atendendo a (41) tem-se agora 


a 
es 1 
Cr - = 2% 2, dy Ami (46) 
m 
Por definição de matriz inversa 
= 
papa 2,0 (47) 
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o que mostra que 
4 = 2) Pim mi (48) 
m 
em concordância com (44), 


É claro que nesta base se compreende que as 
expressões (20) a (23) estão perfeitamente correctas 
desde que se entenda, como está indicado em (22) 
que a transformação B se aplica a v. Considere-se 
na verdade de novo o conjunto de expressões (3) 


A (Vivava) a (Vivava) a, 2,5 2, 
do dog Bog 
dg 2go Rag 
(49) 
= (a, Y% + au ba + Ag Y 224 + Aga La + 


+ Aba AY + Agta T Asa Ya) 


Admitamos que se multiplica o 2.º membro à 
direita pela matriz I(B) não à maneira habitual, o 
que equivaleria a multiplicar a matriz 

(dy Pa 65) 


mas multiplicando as funções “j Obtém-se a matriz 
linha seguinte (que se escreve, por facilidade de es- 
crita, sob a forma de matriz transposta de uma ma- 
triz coluna) 


tr &,; Di Wu + Au Da Ya A, Day Ya + Agy Dia Ya + 
+ az Dao Vo 7 As Dao Va + Ag Dig Lo + Ay Dana + 
+ Ag Das Ya 


2a Dt + Aja Dar Va + Aço Day Vs + Ago ba + 
—+ BagDasba + AsaDasta + Aga Distr + Aga Data + 
+ Saga Dao 


A Di 4 + Aja Day Va + Asa Dyy Ya + Ag Dia by + 
+ AggDoata + AogDasta + Bag Digh + Ag Dasva + 
+ Ass Das ba 


=trgr XD A Fa ZDoy My Ha Da Ra 
Y XD Ag tt ED Aa to Day Aa 


1 | ] 
4 Xby Agra Ebay Ag ts Day Rpg 


= (tata) 2% EP do XPyy dps 
XD da Uboçõno XDoy Ayg 
2% Hr do Soy Rpg 
ou seja, obtém-se, como seria de esperar, 


(Vivava) T(B) DA) 


TÉCNICA 436 


Justifica-se assim que se escreva, agora, e em 
pleno acordo com (23) 


BA (Wubot) = (Wav) POB) MA) (50) 


4 — Podemos agora resumir as conclusões a que 
chegámos. Se considerarmos três transformações A, 
B, C, de um grupo de transformações, podemos de- 
finir as matrizes I(A), T(B) e T(C) associadas a 
essas transformações quando se toma como base um 
dado conjunto de funções v, através de relações de 
forma (1), em que as funções 4 se agrupam forman- 
do uma matriz linha, eventualmente multiplicada à 
direita pela matriz associada a cada uma das trans- 
formações. 

Se entre as transformações A, Be €C existe a 
relação (cfr. (8) ) 


C = BA (51) 


então a matriz associada a €, I(C), que provoca nas 
funções v o efeito equivalente ao de aplicar A às 
funções | e em seguida ao resultado aplicar B, está 
relacionada com as matrizes T(A) e I(B) associadas 
a A e B, respectivamente, e calculadas tomando 
como base as funções |, através da relação 


P(C) = T(A) I(B) (52) 


o que corresponde a multiplicar as matrizes pela 
ordem inversa da aplicação das transformações. 

Uma análise mais profunda do problema permite 
mostrar no entanto que uma vez que se aplica A ao 
conjunto de funções v se obtém um novo conjunto 
de funções &, combinação linear das funções iniciais. 
A matriz associada à transformação B quando a 
base é o cconjunto & é dada por 


P(B) = TI(A) I(B) I(A) (53) 


em que T-1(A) é o inverso da matriz I(A). Nesta 
ordem de ideias a aplicação sucessiva da tranfor- 
mação A e da transformação B às funções Y é tradu- 
zida pela matriz I["(C) dada por 


P(C) = P(A) P(B) = PB) PA) (54) 


sendo agora o produto das matrizes feito pela mes- 
ma ordem pela qual se aplicam as transformações. 
E de salientar que em rigor isto significa, em última 
análise que a transformação B se aplica às funções (, 
iniciais, e não ao resultado da aplicação de A ao 
conjunto «, 


5-—E fácil verificar para um caso concreto que 
na verdade as presentes considerações estão correc- 
tas, Consideraremos por exemplo o caso do grupo C,, 
e as seguintes transformações 


Pol —» rotação de 2X 120º no sentido contrá- 
rio ao dos ponteiros do relógio 

Oy —> reflexão no plano vertical que passa 
pelo eixo dos xx 
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dy —> reflexão no plano vertical que faz um 
ângulo de +120º com o plano corres- 
pondente a o, 

«+ —> reflexão no plano vertical que faz um 
ângulo de 120º com o plano corres- 
pondente a o, 


Se identificarmos 


A com C, 


B com o,, 


da tabela de multiplicação do grupo €,, é possível 
extrair que 


C= BA = ow (55) 
sendo de notar que 
AB = o; (56) 


Por outro lado admitamos que 


4 =x er 
Yy =y er (57) 
Y%y =Z er 


Nessa circunstância, o conjunto de funções dr 
Ps by» definidas em (14) e (29) é 


“= sao SR v3 o a vu Ra 
is va vo 
C. Ya Rd ud RE É er =z- 2 WD va 


Cs 4-2 ecr=u, 


Então, e de acordo com a definição (1) 


na=ncy=/—1/2 —v 8/2 ON 
v 3/2 — 1/2 0 (59) 
0 0 1 


Pode ainda verificar-se que 


MB)=Mo)= / & 
O —i 0 (6) 
& A má 
How)= /—1/2 v32 O 
v 3/2 1/2 O (61) 
WU 0 1 
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P(ory) =  — 1/2 — 3/2 0 
—v 3/2 12 0 (62) 
O ) 0 
MC,)= ac 1/2 v 3/2 0 
—v 3/2 — 1/2 0 (63) 
0) 0 1 
É fácil ver que 
P-(A)j=MC)=  —1/2 —v3 O 
VS] —l/z O (64) 
o o 1 


Atendendo a estes valores para as matrizes po- 
dem agora verificar-se as relações anteriormente 
deduzidas. Assim tem-se 


HA MB= /.4/2 vs O 1 0 0N= 
v3/2 1/2 0 0-1 0 
0 0 1 O 01 
= /—1/2 v3/2 0N =NMC) 
v3/2 1/2 0 cet 
0 0 1 


de acordo com (52). 

Por outro lado a matriz de transformação I'(B) 
de transformação B quando se toma como base as 
funções & é, de acordo com (53) 


MB)=T (A)MBIMA)S / —1/2 V32 0N 
—v 3/2 —1/2 0 
0) 0) 1 
1 oo cc 1iD «cn BIB ONT 
O 1 Q v 3/2 — 1/2 0 
E) 01 0) O z (66) 
É 1/8 050% [==1/3 = 6X 
—v 3/2 1/2 0 v 3/2 1/8 O 
0) 0 1 E) O 1 
o AS". v 3/2 0 
v 3/2 1/2 0 
0) 0 1 


Então o resultado global de aplicação A ao con- 
junto | e depois aplicar B, não ao resultado mas às 
funções | na sua nova configuração é 
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MAB) = 


—1/2 —v8/2 0 
VS - «1/2 0 
0 E 
(67) 
'—1/2 V8/2 0N= —1/2 —vy3/2 0) 
V3/3 1/3 0 DOE 7; ao, 7, 
0 0.4 0 54:32 


Este resultado pode aliás obter-se directamen- 
te de 


O —t 0 
o “027 
(68) 
DE Dani 8/8 ONT 1/2 «8/8. 0: 
V3 —1t/3 0 —v3/2 1/2 0 
0 0 1 , 0 0 1 


o que está de acordo com (54), 

Como é óbvio poderiam arranjar-se outros exem- 
plos, dentro do grupo C€C, ou com outros grupos de 
transformações. 


6 — E curioso observar que estas conclusões a 
que se chegou dependem da forma como é definida 
a matriz associada a uma transformação. Na reali- 
dade, a aplicação de uma dada transformação A de 
um grupo de transformações a um conjunto de fun- 
ções Y que eventualmente irá constituir uma base 
para a representação não define por si só a matriz 
de transformação. O que é essencial garantir para 
que o conjunto | constitua uma base adequada é 
que por aplicação da transformação A a cada uma 
das funções | se obtenha uma combinação linear das 
funções 4. Isto é, deverá ser, se por exemplo, como 
temos vindo a considerar, forem três as funções 1, 


Eos DA 

At = My + dada + Aa Pg 

E Ps: SÃO o 

Av = MW + hs va + Ag %g 

o o o 

Aw = Mw + ha Pa Tr dava 
Estas expressões são perfeitamente equivalentes 
a (2). Simplesmente usou-se a forma (2) porque se 
admitiu que a matriz associada a A, T(A) devia ser 
definida de acordo com (1), isto é, que 1,  %W 
formavam uma matriz linha e portanto I(A) teria 


que aparecer como uma matriz quadrada de ordem 3, 
colocada à direita de 


(ri va Va) 


Nessas circunstâncias, sendo a,, o elemento ge- 
nérico da matriz T(A) e usando a convenção habi- 
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tual de o índice (k) representar a coluna e o índi- 
ce (i) representar a linha dessa matriz a forma (1) 
de definição de T(A), 


Aw= Sa, as 
k 
(70) 
implica a identidade de (2) e (69) ou seja 
D 
A; = 4 
(71) 


Daqui decorrem todas as expressões anteriormen- 
te consideradas. Em particular aquelas que permitem 
obter o conjunto & a partir do conjunto | e onde 
se manteve a convenção de considerar v,, “, Y% e logo 
Pp Po 9) constituindo matrizes linha (cfr. (35) 
e (36) ). 

Suponhamos que se considera agora que “4, “., 0: 
se agrupam formando uma matriz coluna. Nessa al- 
tura a matriz associada a A deverá ser igualmente 
uma matriz quadrada, de ordem 3, mas que multi- 
plica à esquerda a matriz coluna 


! (72) 


É fora de dúvida que esta nova matriz deva es- 
tar relacionada com a matriz P(A) já definida ante- 
riormente. Para vermos como as matrizes associadas 
a A de acordo com as duas possibilidades de agrupar 
vi Y% € & estão relacionadas representemos por a, 
o elemento genérico da nova matriz. Ter-se-á agora, 
em vez de (1) (6) 


Bad x 4 
Avy= Sant 


(73) 


Desenvolvendo para os vários valores de y, tem- 
-se agora, em vez de (2) 


Ad = ay + ão + ada 
Ava = et + Cato + Zagoy (13) 


| e | J ' 
Av, = Gyvi + Kara + Zura 


Em termos das relações (69) tem-se agora 


(75) 


MO) = amu 


o que permite desde já escrever, comparando com 
(71), que 


(76) 


da = aim 
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ou seja que à nova matriz associada à transforma- 
ção A é a matriz transposta de T(A), isto é, a 
matriz que se obtém trocando as linhas e colunas 
de T(A). 

Esta conclusão é aliás evidente se se atender a 
que a matriz transposta de um produto é o produto 
das matrizes transpostas escritas em ordem inversa, 
e que a matriz transposta de uma matriz linha é uma 
matriz coluna, 

Na realidade, se aceitarmos a definição (73) de- 
verá escrever-se, em vez de (3) 


/ o —— 4 - 
A ra = &j Yo Gig Y ns; 
y 4 
2 Go Log Gas é 
' 
Ya Ga Tg Rg |) Ya 


- 2% + ãxata + mata 
Zu TF Mata + ab (17) 
ant + Mova + mas ba 


Se usarmos as matrizes transpostas das indica- 
das em (3) deverá ter-se 


MEAN ds & did 
Atr(WvVa) = tr | (MV) a, dd, dd 


dd; do das 
(78) 
ou seja 
A UN Etr a, dd ds x / W = 
Ya do do dos Yo 
f 
na da dy Ass Ya 
a A MM ds | Ya 
| 
do do dos Ya (79) 
2, ds, Egg / 4 
o que mostra, por identificação com (77), que 
&s Gs Gis = tr d, do &, 
Qo) oz os do) Ao do; (80) 


&a ga Gyg a a 


em perfeito acordo com (76). Seguindo a convenção 
habitual usaremos 


tr MA)J=IMA) (81) 


E claro que para duas outras transformações 
Be Cas matrizes associadas já não serão as matri- 
zes T(B) e T(C) definidas em (4) e (5) mas sim 


as respectivas matrizes transpostas T(B) e TTO). 

Para saber como estas matrizes estão relaciona- 
das quando entre as transformações A, Be C existe 
a relação (8) (ou (51) ) podem repetir-se os racio- 
cínios anteriormente feitos. A resolução do problema 
simplifica-se no entanto se usarmos as matrizes trans- 
postas em (52) 


ITÉ) = 17B) ITA) (82) 


Este resultado mostra que quando se usa como 
matriz associada a uma transformação aquilo que é 
definido pela relação (73) e se consideram três trans- 
formações A, Be C entre as quais existe a relação 


C = BA (83) 


então a matriz associada a C, fio), que provoca 
nas funções y o efeito equivalente ao de aplicar A 
às funções i e em seguida ao resultado aplicar B 


estã relacionada com as matrizes STA) e ÍTB) 
associadas a A e B, respectivamente, e calculadas 
tomando como base as funções v, através da relação 
(82) o que corresponde a multiplicar as matrizes 
pela mesma ordem pela qual as transformações são 
aplicadas. 

Como é evidente poderia chegar-se a esta con- 
clusão repetindo as considerações traduzidas pelas 
equações (9) a (13) tomando como base a definição 
(73) em vez de (1) ou ainda usando as expressões 
homólogas das expressões (14) a (19). 


7 — Mas tal como foi frisado anteriormente esta 
análise não toma em consideração o facto de que por 
aplicação de A ao conjunto y se obter um conjunto q 
que ao ser tomado como base para definir a matriz de 
transformação associada a B conduz a uma matriz 
diferente daquela que se obtém quando a base é o 
conjunto y. Agora entre o conjunto & e o conjunto 1 
existe o conjunto de relações (74). Ou seja, sob forma 
matricial deve escrever-se, em vez de (35) e (36) 


P, — A) vi 
Po Ya (84) 
Pg Ya 
e 
7 mm a y 
| = [ra] 
“a Po 
va Ps 
= rã (A) Pay 
Pa (85) 
Pa 
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Estas expressões podem obter-se directamente de 
(35) e (36) transformando as matrizes correspon- 
dentes. Do mesmo modo se obtêm as expressões ho- 
mólogas de (38) e (39) 


ai Adi, = ITA) NB) 
da wu |(86) 
Ps da 
P/ AN =yrmbema/ * 
A A Zé, 
Ps Ps 


resultado que mostra que a matriz associada à trans- 
formação B quando as funções base são as funções q 
e se usa a definição (73) é 


md mad mad Ed 
P(B) = IA) IB) T-I(A) (88) 


Este valor da matriz associada a B conduz à 
expressão homóloga de (44) 


TC) = ITA) ITB) (89) 


quando se substitui P(B) por rTB) em (82). 

Tal como sucedia no entanto com a dedução tra- 
duzida pelas expressões (20) a (24) este resultado 
só é válido quando se admite que a transformação B 
actua sobre o conjunto v, É instrutivo repetir aqui 
esses raciocínios. 

Seja com efeito q; o elemento genérico da ma- 
triz associada a A obtida com base na definição (73) 


As! : Ta y 
4 > ik (90) 


Para a transformação B será igualmente, sendo 
8; º elemento genérico 


E q (91) 


(92) 
ou, atendendo a (91) 
BAS Wa WB 4% 
pa 
i K dr NS 3 
= (Dx Pay) 4 
= ik Cj) 73 
(93) 
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Por definição das matrizes P(A) e ÉTB) tem-se 
então 


Ec A = Ta) FB) á 
Yo Yo (94) 
N va Ya 
ou seja 


TO) = ITÃ) ITB) 
(95) 
como em (89). 


8 —Para aplicar estas expressões a um caso 
concreto pode usar-se como anteriormente, o caso 
do grupo €, e o conjunto de funções v v, vw defi- 
nido em (57). As relações (58) e as relações homó- 
logas que se podem obter para as outras transforma- 
ções continuam evidentemente a ser válidas, Simples- 
mente as matrizes associadas às várias transforma- 
ções serão as transpostas das anteriormente apresen- 
tadas, como se viu com toda a generalidade em (80). 


Seja então 


— 1/2 +v 3/2 0 


NÃ) = PC?) 


= 88  —=1/2 0 
0 0 1 
(96) 
1 0 0 
P(B) = Pow) = 
O —1 0 (97) 
0 0 1 
Tem-se agora 
TO) =TTB) ITA) = 
o UR 
Oo 01 
(98) 
= —l/2 3/2 0) = -—t/B 8/8 O 
87 t/B -QO v3/2 1/2 0 
0 0 3; 0 0 1 


de acordo com 


ovv CS = ow 
(99) 
como já indicado em (55). 


9 —Como se vê, quando a matriz associada a 
uma dada transformação é definida por uma expres- 


são da forma (73) a matriz que traduz o efeito da 
aplicação sucessiva das transformações A e B a um 
conjunto de funções | é o produto das matrizes cor- 
respondentes a A e B pela mesma ordem em que as 
transformações se aplicam. Essa ordem é invertida 
quando as matrizes associadas às transformações A, 
B e C são obtidas através de expressões da forma (1). 

Assim sendo, é o conjunto de matrizes que se 
obtém com base na definição (73) que apresenta uma 
tabela de multiplicação idêntica à do grupo de trans- 
formações a que A, B e € pertencem. É portanto a 
definição (73) que deve ser escolhida quando se pre- 
tende determinar o grupo de matrizes isomorfo do 
grupo de transformações. 

Pode argumentar-se aqui, na linha do que foi 
apresentado anteriormente, que uma vez que a apli- 
cação da transformação B incide sobre um conjunto 
de funções que não é o conjunto inicial o que se 


deve usar é, não a matriz P(B) (ou 1B) ) mas sim 


a matriz I'(B) (ou TB) ). A ser assim, como se 
viu, deveria ser o conjunto de matrizes obtido com 
base na definição (1) que deveria ser usado, já que é 
esse que conduz a um produto das matrizes I(A) 
e T(B) pela ordem com que as transformações são 
aplicadas (cfr. (54) ). Este argumento é contudo fa- 
lacioso. Em primeiro lugar porque o problema que 
está em discussão é o de saber qual dos conjuntos 
de matrizes que se obtém tomando como base as 


funções y (usando a definição (1) ou a definição 
(73) ) constitui um grupo isomorfo do conjunto de 
transformações em estudo. Em segundo lugar é de 
salientar que se se aceitasse o argumento seria ne- 
cessário dispôr de uma matriz representando a trans- 
formação B para cada transformação A que se apli- 
casse inicialmente (cfr. (40) e (88) ). E a aplicação 
de nova transformação exigiria novo conjunto de 
matrizes. 

Se bem que os raciocínios que conduzem às equa- 
ções (40) e (88) estejam perfeitamente correctos 
elas não respondem ao problema proposto. É na ver- 
dade o conjunto de matrizes obtidas com base em 
(73) que constitui o grupo isomorfo pretendido, como 
aliás se verifica quando se considera um caso con- 
creto (cfr. eg. (96) — (99) ). 
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Transmissão de ruído através dum dispositivo 
retentor de amostras !) 


RESUMO 


Com base num modelo apropriado dum disposi- 
tivo retentor de amostras, estabelece-se a expressão 
da função densidade espectral à saída, devida à pre- 
sença à entrada, de flutuações com função de auto- 
-correlação especificada, 

O modelo utilizado para representar o retentor 
leva em linha de conta os valores não desprezáveis 
da constante de tempo e do tempo de captação do 
dispositivo. 

Apresentam-se os resultados de medidas efectua- 
das num dispositivo real, para confirmação das pre- 
visões teóricas e ilustra-se a aplicação prática do re- 
tentor, com os resultados obtidos pelo seu uso na re- 
cuperação e reprodução da forma de curva duma 
grandeza variável imersa em ruído. 


1. INTRODUÇÃO 


O retentor de amostras (1) é um elemento de 
memória do domínio analógico, com larga aplicação 
na técnica do comando automático discreto [1] (2) 
e nas medidas eléctricas em geral. 

Neste último domínio, encontram-se os diversos 
métodos de registo gráfico de fenómenos periódicos 
ou repetitivos baseados na amostragem, permitindo 
a reprodução da sua forma como funções do tempo, 
utilizando instrumentos registadores com reduzida 
largura de banda, sejam eles oscilógrafos de raios 
catódicos (3) ou até registadores potenciométricos, [3]. 

Na referência [3], descreve-se um método de 
identificação de sistemas baseado precisamente na 
amostragem, cujo desenvolvimento motivou o pre- 
sente estudo. 


J. F. BORGES DA SILVA 
Doutor-Engenheiro 1.S.T. 

Prof. Extraordinário do 1.S.T. 
J. M. VIEGAS GONÇALVES 
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ABSTRACT 


Based on an appropriate model of a sample-and- 
“hold device, a theory is developed leading to an ex- 
pression for the output spectral density function due 
to the presence at the input of fluctuations with a 
specified auto-correlation function. 

The sample-and-hold model used takes into 
account the non-negligible time constant and acqui- 
sition time of the device. 

Results of measurements made on a real device 
are presented to confirm the theoretical predictions 
and an applications example is described being illus- 
trated with results obtained, using a sample-and-hold 
device, for the recovery and display of wave-forms 
buried in noise. 


A finalizar, cite-se ainda o emprego do retentor 
na obtenção de valores instantâneos de grandezas 
variáveis sob forma digital para transmissão directa 
a um computador, em sistemas de medida com aqui- 
sição e processamento automático de dados [4]. 

Nas aplicações referidas acima em que o reten- 
tor é utilizado para a reprodução de formas de curva 
com forte redução de largura de banda, surge como 
propriedade de grande interesse a possibilidade de 
imunidade em face de ruído. Com efeito, neste género 
de aplicações, o retentor comporta-se em relação à 
grandeza periódica a reproduzir, praticamente como 
um detector síncrono. Filtrando a saída, e sem afec- 
tar a reprodução da forma de curva que exige apenas 
uma pequena largura de banda, é possível rejeitar 
os efeitos de flutuações não correlacionadas sobre- 
postas na entrada, permitindo a reprodução fiel, mes- 


() Publicação do Laboratório de Medidas Eléctricas do I. S. T. integrada no Projecto de Investigação TLE-5 subsi- 


diado pelo Instituto de Alta Cultura, 
() Sample-and-hold device. 
() Sampled data control systems, 
() Sampling oscilloscopes. [2]. 
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mo no caso de grandezas completamente submersas 
em ruído. 

O grau de rejeição possível, como é natural, 
varia inversamente com a largura de banda do filtro 
empregado e é portanto proporcional ao tempo dis- 
ponível para efectuar a reprodução, 

Sendo o retentor fundamentalmente um interrup- 
tor comandado, é um dispositivo linear. As respostas 
à grandeza a reproduzir e ao ruído adicionam-se sim- 
plesmente. Nestas condições, o estudo da transmissão 
unicamente de ruído através do retentor é suficiente 
para caracterizar este aspecto do seu funcionamento. 

O objectivo central do presente estudo consiste 
assim em partir dum modelo realista do retentor com 
consideração dos valores finitos da sua constante de 
tempo e tempo de captação, e com base nele prever 
por via teórica a densidade espectral das flutuações 
à saída quando excitado à entrada por flutuações 
estocásticas com função de auto-correlação bem de- 
finida, submetendo os resultados à confirmação ex- 
perimental. Ilustra-se ainda uma possível aplicação 
prática deste estudo com os resultados obtidos na 
recuperação e reprodução de formas de curva de gran- 
dezas imersas em ruído usando um retentor. 


2. DESCRIÇÃO DO MODELO DO RETENTOR 


O modelo adoptado está esquematizado na Fig. 1, 
v consiste essencialmente num interruptor comanda- 
do S seguido dum circuito integrador RC. A entrada 
é representada pela tensão y(t) e a saída pela ten- 


yít) Cc x(t) 


To= Rolo - 


Fig. 1 
Modelo do retentor 
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são x(t) aos terminais do condensador, supostos em 
vazio. Com o interruptor S aberto a tensão de saída 
x(t) permanece constante com o valor que tinha no 
instante de abertura. 

O interruptor S é fechado a intervalos regula- 
res T (período de amostragem) e permanece fechado 
durante um breve intervalo de duração T. (tempo 
de captação). Durante este intervalo o circuito RC 
«adapta-se» ao valor da tensão de entrada y(t) im- 
posta, em maior ou menor grau consoante a sua 
constante de tempo T,=R, C, permitir. 

Idealmente deve ser T.<T e T,<T, a fim da 
tensão x(t) retida no condensador representar prati- 
camente um valor instantâneo de y(t). Na teoria que 
se expõe no seguimento, admite-se ser verdadeira a 
primeira condição, mas não necessariamente a segun- 
da, sendo os valores finitos de T, e T, levados em 
consideração. 


3. ANÁLISE DO PROCESSO DE CAPTAÇÃO DUMA 
AMOSTRA 


Com o interruptor S da Fig. 1 fechado, tem-se 
em valores instantâneos, 


dx 
ro Pd ada (3.1) 
onde T, =R,C, representa a constante de tempo do 
retentor. 

Seja f(t) um integral particular de (3.1) corres- 
pondente à tensão y(t) imposta e h(t) o integral ge- 
ral de (3.1) com o segundo membro nulo, Nestes ter- 
mos a solução de (3.1) é da forma, 


t—t 
X(O=f(0) +HhCO=f(0+h, exp | - | 
To 


Tomando para t o instante em que se inicia a 
captação da amostra n e designando por x. o valor 
retido no final da amostragem anterior, Fig. 2, tem- 
-se, calculando o valor da constante h que assegura 


a continuidade de x(t) em t=t, e para t>t, 


É — t 
x(t) = f(t) -[4 =X,. leso [- | 
l Es 


Designando por x o valor de x(t) retido no final 
da amostragem n, tem-se, fazendo t=L+T, na ex- 
pressão anterior, 


X, = f(t, + 2.4) nisto [H(t) | ET (3.2) 


pe 
Op = €Xxp - | 


A expressão (3.2) permite calcular a sucessão 
de valores x das amostras retidas a partir de su- 


onde, 


(3.3) 
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cessivos valores de f(t). Note-se que esta função não 
é mais que a resposta do circuito R,C, da Fig. 1 à 
tensão y(t) imposta, com o interruptor S permanen- 
temente fechado. 


f(t) 


Fig. £ 
Pormenor da captação duma amostra 


4. DENSIDADE ESPECTRAL RESULTANTE DA AMOS- 
TRAGEM DE FLUTUAÇÕES ESTOCÁSTICAS 


Admitir-se-á que o intervalo de captação T é 
sempre muito pequeno comparado com o período de 
amostragem T, Fig. 1. 

LT (4,1) 

Nestas condições, a tensão x(t) na saída do re- 
tentor pode ser considerada como uma sucessão de 
impulsos rectangulares justapostos de duração T, 
Fig. 3, cujas amplitudes são dadas pelos valores da 
sucessão x. Em conformidade escreve-se, 


+50 
x(t) = >, x, g(t—nT) (4.2) 
— DO 
onde, 
=L,Ogr<T 
4.3 
ad DE tr (4,3) 


representa um impulso de duração T, e amplitude 
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unitária, cuja transformada de Fourier é dada por, 


T=T 
ro [5 g(r) exp (— Jor) dr=' 


exp(— jwT/2) (4.4) 


(o T/%) 

Para efeitos desta análise, a tensão y(t) imposta 
é encarada como uma variável estocástica, bem como 
a sucessão de valores x resultantes da amostragem 
e consequentemente também a tensão x(t) à saída 
definida por (4.2), 

A análise, pormenorizada no Apêndice-l, conduz 
a uma expressão para o valor expectável da densi- 
dade espectral de x(t) (por Hz de largura de banda), 
da forma, 


o = 


Co 
- 
P(o)= | lo) | Xo +2 2 X, cos(suT) (4.5) 
B=: 
onde 


x. = <XX > (4.6) 

representa o valor expectável da auto-correlação da 

sucessão x, para um atraso de s=1,2,3,... amostras. 
Numa vizinhança da frequência zero tal que, 


o € (27 /T) (4.7) 


pode considerar-se a densidade espectral praticamen- 
te constante com o valor igual a P(0), Tomando o 
limite (1) -> 0 em (4.4) obtém-se de (4.5), 


So 
P(0)=2T|X+? ps x. (4.8) 


resultado a aplicar no cálculo do valor quadrático ex- 
pectável das flutuações à saída dum filtro de banda 
estreita colocado após o retentor. Designando essa 
saída por z(t) será, 


tn.a tn.1 tn 


tnes tnea t 


Fig. 3 
Saída do retentor 


2 
L=0,P(0) (4.8) 


onde s), designa a largura de banda eficaz para ruído 
branco do filtro, expressa em Hz. 


5. EQUAÇÃO DE RECORRÊNCIA DA SUCESSÃO X. 


Tome-se a equação de recorrência da sucessão 
x, (3.2) escrita na forma, 


X— do = TT) — go ft) 


referida à amostra n cuja captação se inicia em t=t,., 
e ainda a mesma equação referida à amostra n-s 
cuja captação se inícia em t=t -sT, 


x — 9 X 
ns O nes-l 


=f(t—sT+4+T)-d fit —sT) 
n ' 0 um 


Multiplicando membro a membro estas equações 
e extraindo o valor médio, obtém-se em termos de 
valores expectáveis a equação, 


Ke — (do + 6, )X, Fhega =P = 14) = 


— (d+ 6,) FIST) + F(ST HT) (5.1) 


onde a função F'(r) designa o valor expectável da 
auto-correlação de f(t), isto é, 


F(T) = <f(t) flt-m)> 


A equação (5.1) é uma equação de recorrência 
de segunda ordem em X,, linear de coeficientes cons- 
tantes que permitirá obter a sucessão de valores X. 
pretendidos, dada a função F'(r) de auto-correlação 
de f(t). Esta função, como se salientou no final de 3,, 
não é mais do que y(t) filtrada por acção da cons- 
tante de tempo T, = R,€C, do retentor. 


6. SOLUÇÃO DA EQUAÇÃO DE RECORRÊNCIA 


Para exemplificar o método de cálculo, consi- 
dera-se o caso elementar de y(t) ser ruído branco 
filtrado por acção dum simples filtro RC passa-baixo 
de constante de tempo T.. 

O caso mais geral dum filtro com função de 
transferência racional pode ser tratado por simples 
extensão deste mesmo método. 

Nas condições enunciadas, e desde que seja, 


TAT, (6.1) 


obtém-se. (Veja-se Apêndice-2), 


E(=) = x exp(- |z|/To) +a, exp(— IxI/T,) (6.2) 


onde, 
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à À 2 
ca DP na Es: 6.3 
04 Ty or sp 2 ? em 
0 1 1 
Na escrita da equação de recorrência (5.1) para 
este caso particular introduza-se ainda, além de q, 


já definido em (3.3), 


d = exp( — TT) 
(6.4) 
A = exp( —'T/T) 


considerando-se em primeiro lugar as equações rela- 
tivas a |S| > 1, em que os elementos da sucessão X 
intervenientes, ou são todos de índice positivo, ou 
todos de índice negativo. 

A substituição de (6.2) na equação de recorrên- 
cia (5.1) conduz neste caso a, 


SDL Xa—(do+ IX AX =Ball (65) 
onde, 


by = [ (64 ) cod 6) la (6.6) 


Admitindo que, 


da E A, 


o1 seja, atendendo às definições (3.3) e (6.4), 
T/T TT, (6.7) 


a equação (6.5) admite a solução particular da forma, 


AA” 
» háso | 
onde a constante A, se pode obter por substituição 


em (6.5), e indiferentemente para s>1 ou para 
s< —l, com o resultado, 


(6,48) — (do + 87) 


À, = 14 (6.9) 
(+27) — (do +87) 

A esta solução particular pode sempre sobrepor- 

-se uma solução da equação homogénea correspon- 

dente a (6.5) com o segundo membro nulo, São possí- 

veis dois tipos de solução independentes da forma, 


s “Ss 
*, 'o 
Atendendo a que a sucessão X, representa uma 
auto-correlação, conclui-se que deve ser função par 
do índice s e tender para zero com |s| —» c0. Como 
bo <1, a solução q s aplicar-se-á no domínio s <1 e 
a solução &,-s no domínio s < -1, pelo que a solução 
geral de (6.5) deve ser tomada na forma. 


£,=B a” vas (6.10) 


TÉCNICA 436 


onde B é uma constante arbitrária. O valor desta 
constante deverá ser fixado de modo a que esta solu- 
ção satisfaça ainda a equação de recorrência (5.1) 
para s = 0. Nesta equação surgem simultaneamente 
X. e X, cujos valores devem ser iguais, pelo que a 
equação toma a forma, 


2X, — (do + 8, )Xo =2E(TÇ) — (go +87 )E(M) 


onde, por substituição de (6.10), se obtém o valor da 
constante, 


1—| (ho +48! 2) — 
fo 


B = «q + 


(6.11) 
— (do + 6 — 2)A, 


O resultado expresso por (6.10) mostra que a 
correlação dos valores das sucessivas amostras tem 
duas origens distintas. Por um lado, como era fácil 
de prever, resulta dessa correlação existir já nas 
flutuações impostas na entrada do retentor. Este 
efeito corresponde à parcela em A, de (6.10), 
que será apreciável desde que a constante T que 
caracteriza o filtro de entrada seja comparável ao 
período de amostragem 'T, veja-se (6.4). Por outro 
lado, a constante de tempo T, do retentor, se for 
comparável ao tempo de captação T,, determina 
igualmente a existência de correlação das sucessivas 
amostras através da parcela em gs de (6.10), 
veja-se (3.3). O facto do regime livre associado à 
constante de tempo T, = R,C, do retentor não ter 
tempo de se extinguir totalmente durante o intervalo 
de captação, faz com que, independentemente dou- 
tras circunstâncias, o valor de x, da amostra n de- 
penda do valor x. da amostra anterior (3.2), con- 
tribuindo para a correlação das sucessivas amostras, 


7. CÁLCULO DA DENSIDADE ESPECTRAL A SAÍDA 


O facto de X, ser dado por uma soma de pro- 
gressões geométricas de razões g, e À, (6.10) per- 
mite somar a série (4.5). Representando a função 
cos(SwT) na forma exponencial, 


cos(swT) = 5 | expliso?) + exp(—jsoT) | 


obtém-se (4.5) como uma soma de progressões geo- 
métricas de razões complexas dadas pelas exponen- 
ciais de expoentes, 
Yo =—(T/T) +iwvT 


LE T/T) +joT (71) 


e os seus conjugados. 
Fazendo uso da expressão (4.4) para T(y) O 
resultado da soma exprime-se por, 
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sen?(wT/2) 


P(v) = 2T (9/2) 


Rel B coth(-y9/2) + 


(7.2) 
+ A coth(-1/2) | 


Em particular, para w—>0, 


P(0) =2T |B coth(T, /2To) + A coth(T/2T,) | (7.3) 


onde as constantes A e B são calculadas segundo 
(6.9) e (6.11) em função das razões de tempos 
T/To T/T, e T/T,. 


8. RESULTADOS NUMÉRICOS 


Para obter uma ideia do andamento da função 
descrita por (7.2), considere-se o caso particular dum 
retentor ideal em que T, —> 0. Nestas circunstâncias 
desaparece a parcela de coeficiente Be A >Y2, 
obtendo-se, 


sen?(w 1/2) 
(wT/2)? 

Se além disso admitirmos que T <T, ou seja, 
que o tempo característico da correlação à entrada 
é desprezável face ao período de amostragem, o fac- 
tor em coth( —y, /2) fica praticamente constante e 
igual à unidade. Nesta situação limite o andamento 
da função P(«w) é traduzido pela função (senx/x)* 
que é assim característico das situações em que as 
amplitudes das sucessivas amostras não apresentam 
correlação apreciável. Podem então encarar-se os 
efeitos da correlação como uma modulação deste es- 
pectro básico por efeito dos termos da expressão 
(7.2) dependentes de y, e y,. 

Para ilustrar o seu efeito considere-se o termo 
dependente de y,, 


P(9) = 2T R e coth(-7,/2) | Cs (8.1) 


1l-—à 
1 


Re 4 coth(-y,/2) ! E A RR SE 
| 1+1 — 22,cos oT 
1 


Quando multiplicado por um factor positivo como 
em (8.1), os máximos ocorrem nos zeros do primeiro 
factor e os mínimos a meia distância entre os má- 
ximos, provocando um achatamento das arcadas la- 
terais da função |senx/x| e reforçando a arcada 
central, Fig. 4. 

Os resultados que se obtém para várias combi- 
nações dos tempos característicos encontram-se re- 
produzidos nas Figs. 5 e 6. 

Nestes gráficos normalizou-se a escala das fre- 
quências tomando a frequência de amostragem 
f =1/T para referência, fazendo 


wT =2 q(f/£,) 


e usando o cociente adimencional f/£, em abcissas. 
Quanto às ordenadas usou-se a quantidade adi- 
mensional definida por, 
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1 
FE 


p= (P/YaT) =(PE) [Yo 


6 | mé 


Esta quantidade é proporcional ao valor eficaz das 
flutuações à saída e não ao valor quadrático das 
mesmas, sendo o valor eficaz das flutuações à en- 
trada tomado como referência. Corresponde assim 
mais directamente aos valores que efectivamente são 
medidos nos ensaios experimentais, A introdução do 
factor 1/T=f, refere o resultado à largura de banda 
de f Hz em vez de 1 Hz com o que se obtém um 
resultado independente do valor absoluto do período 
de amostragem T, que somente intervem através da 
razão T/T.. 


(tl 


p 
Ver 


Fig. 4 
Deformação do espectro básico. (a) v2!sen x/x |, x=wyT/2. 


(b) [Re (coth (—y,/2))] e? 


No primeiro grupo de resultados, Fig. 5, é T,= 
=0,1T.. A constante de tempo do retentor pode con- 
siderar-se curta, a correlação das amostras é devida 
praticamente apenas à coerência da tensão de entrada 
expressa pelo tempo T| e é importante na medida 
em que este tempo for comparável ao período de 
amostragem 'T. A curva relativa a T =0,02T traduz 
praticamente a situação limite de ausência de corre- 
lação, apresentando a forma característica da função 
|senx/x|. O progressivo achatamento das arcadas 
com o aumento de T, é evidente. 

No segundo grupo de resultados, Fig, 6, fez-se 
T,=T, de modo a por em evidência os efeitos da 
constante de tempo do retentor. A deformação das 
arcadas é em geral mais acentuada e subsiste mes- 
mo para pequenos valores de T,, isto é, mesmo quan- 
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do a coerência da tensão à entrada deixa de ser im- 
portante. 


9. RESULTADOS EXPERIMENTAIS 


9.1 — Densidade espectral das flutuações à saida 
do retentor. 


Para verificação dos resultados da teoria expos- 
ta, realizou-se a montagem cujo esquema se indica 
na Fig. 7, em que: 


Filtro 1 — Filtro RC passa baixo de constante 
de tempo T, (regulável de 50 a 200 us). 
Filtro 2 — Variante a): Analisador de espectro 
com largura de banda fixa (100 Hz) e fre- 
quência central variável comandada por ten- 
são (50 Hz a 1 MHz). Saída proporcional ao 

valor eficaz. 


Variante b): Filtro RC activo passa-bai- 
xo com característica do tipo Bessel de 2.º 
ordem e frequências de corte de 50, 100 e 
200 Hz. Permite obter a densidade espectral 
na vizinhança da frequência zero, 


Al,A2,A3 — Amplificadores lineares (0 a 1 MHz). 


V1,V2,V3 — Voltímetros quadráticos empregando 
conversores termo-eléctricos (5 Hz a 2 MHz). 
Indicação proporcional ao valor eficaz. 


GR — Gerador de ruído com densidade espectral 
constante até 500 kHz, 


Retentor — Dispositivo construído para a amos- 
tragem na identificação de sistemas pelo mé- 
todo da referência [3]. Constante de tempo 
T,=1,6 us (fixa). Tempo de captação regu- 
lável T,=0,8 a 3,2 us. Período de amostra- 
gem T definido externamente. 


O traçado da função densidade espectral é obtido 
directamente no registador, aplicando a saída do ana- 
lisador em ordenadas e a tensão de comando da fre- 
quência central em abcissas. Variante a). 

Não obstante a pequena largura de banda do 
registador, verificou-se ainda a necessidade de filtra- 
gem passa-baixo da saída do analisador a fim de 
reduzir as flutuações a um nível razoável. Utilizou-se 
um filtro RC de constante 3 s, e realizou-se o varri- 
mento da frequência lentamente com durações de 
cerca de 30 minutos. 


A calibração absoluta das ordenadas fez-se apli- 
cando ao sistema de medida ruído de características 
conhecidas. 

Os resultados, tal qual foram obtidos no regis- 
tador, encontram-se reproduzidos nas figuras 8, 9, 10 
e 11, tendo-se sobreposto nestes gráficos pontos cal- 
culados para cada caso com o auxílio da expressão 
teórica (7.2). Estes exemplos correspondem a quatro 
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1 


Fig. 5 


Funções densidade espectral com constante de tempo do retentor curta 


combinações distintas dos valores das razões carac- 
terísticas T/T, e T,/T, com valores fixos de T, e T, 
conforme se indica nas legendas. 


O acordo dos resultados experimentais é, duma 
forma geral, muito bom excepto na vizinhança dos 
zeros que não se observam nos resultados experimen- 
tais, Tal facto deve-se, por um lado à largura de 
banda finita do analisador, por outro lado a imper- 
feições da comutação no retentor. As transições de 
nível entre as sucessivas amostras não são perfeita- 
mente abruptas como se supõe na teoria surgindo 
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além disso pequenas perturbações transitórias nos 
extremos do intervalo de captação, Fig. 17. Estes 
factos determinam a presença à saída dum espectro 
de riscas espúrio com a frequência de amostragem 
e suas harmónicas que vão precisamente incidir nos 
zeros da função densidade espectral prevista. 

A determinação da densidade espectral na vizi- 
nhança da frequência zero fez-se com a variante b) 
da montagem da Fig. 7, estando os resultados mar- 
cados sobre o gráfico teórico da fig. 12, obtido com 
o auxílio da expressão (7.3). 
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Fig. 6 
Funções densidade espectral com constante de tempo do retentor longa 
9.2 — Aplicação do retentor na recuperação de for- RLC — Circuito a ensaiar, cujo esquema está in- 
mas de curva submersas em ruído. dicado na Fig. 14. 

Como exemplo de aplicação, considerou-se o pro- Unidade de atraso — Permite fixar o instante de 
blema da recuperação da resposta impulsiva dum cir- amostragem em cada período da tensão a 
cuito RLC em presença de ruído de amplitude bas- observar com um atraso r em relação a um 
tante superior. Para o efeito foi utilizada a monta- instante de referência, Fig. 15. Esta unidade 
gem representada na Fig. 13, onde além do material possui uma saída analógica proporcional a 7, 
já referenciado a propósito da montagem da Fig. 7 se permitindo a reprodução da amplitude das 
tem ainda: sucessivas amostras em função de 7, (repro- 

dução da forma de curva em tempo expan- 

GI — Gerador de impulsos para a excitação pe- dido). E possível expandir 0,5 us de tempo 

riódica do circuito RLC. real até 1 s de tempo expandido. 
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pY2 (my Hz!) 


000 VALORES CALCULADOS 


Fig. 9 


10 15 kHz 


Função densidade espectral. Resultado experimental obtido para (T,/T,)=0,5 e (T,/T)=1,0 


Y1 — Registo oscilográfico em tempo real da 
tensão amostrada. 


Y2 — Registo oscilográfico em tempo expandido 
do resultado da amostragem, (base de tem- 
po ur). 


A frequência de repetição do fenómeno a obser- 
var foi fixada em 2 kHz, correspondendo-lhe portanto 
um período de amostragem T=500 us (2000 amostras 
por segundo). No oscilograma (A) da Fig. 16 mos- 
tra-se o sinal periódico sem ruído cujo valor eficaz 
era de 30 mV, representando o oscilograma (C) da 
mesma figura a forma de um período em pormenor 
(tempo real). 

O ruído adicionado consistiu na saída de banda 
larga do gerador de ruído filtrada pelo filtro-1, 
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Fig. 13, de constante de tempo T =100 «us, resultando 
em flutuações à saída do filtro com o valor eficaz 
Y,,=200 mV. A relação sinal/ruído à entrada do 
retentor era assim de -16 db. 

O oscilograma (B) da Fig. 16 mostra o aspecto, 
obtido com um único varrimento, da sobreposição do 
sinal do oscilograma (A) com o ruído, tal como se 
apresentam à entrada do retentor. O oscilograma (D) 
da Fig, 16 mostra a forma da curva recuperada obti- 
da por filtragem da saída do retentor pelo filtro-2 de 
largura de banda eficaz 2,=0,4 Hz e em função da 
base de tempo expandida ur. O oscilograma da Fig. 17 
representa um segmento, obtido com um único varri- 
mento, da tensão à saída do retentor antes da filtra- 
gem, (comparar com a Fig. 3). 

Utilizando os resultados do gráfico da Fig. 12, 
obtido a partir da expressão (4.8) e normalizado 
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Fig. 7 
Montagem para o traçado experimental da densidade espectral 


pY? (my Hz) 


DO O VALORES CALCULADOS 


0 5 10 15 kHz 
Fig. 8 
Função densidade espectral, Resultado experimental obtido para (T/T)=0,6 e (T,/T)=0,25 
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00 O VALORES CALCULADOS 


10 15 kHz 


Fig. 10 
Função densidade espectral. Resultado experimental obtido para (Ty/T,)=2,0 e (T,/T)=0,25 


como se explicou em 8., obtém-se para relação entre 
os quadrados dos valores eficazes do ruído à entrada 
Y., e à saída filtrada do retentor Z,,, 


Z ivº =oT po 
ur | ef aa Pi ) 


Atendendo às condições da experiência obtém-se 
uma relação de 4x10-+ ou seja de —34 db, o que signi- 
fica uma relação sinal/ruído de +18 db no oscilo- 
grama (D) da Fig. 16. 

Note-se que os 250 as de tempo real represen- 
tados levaram cerca de 50 s tempo expandido a obter, 
pelo que a recuperação da forma de curva correspon- 
deu ao processamento de cerca de 100 000 amostras. 


10. CONCLUSÕES 


Estabeleceu-se um método de cálculo da densi- 
dade espectral das flutuações à saída dum dispositivo 
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retentor de amostras, partindo da especificação da 
função de auto-correlação das flutuações à entrada 
e com base num modelo do dispositivo consistindo 
num interruptor seguido dum circuito integrador RC. 

Revelaram-se como parâmetros relevantes do dis- 
positivo o período de amostragem, o tempo de capta- 
ção e a constante de tempo do circuito integrador, 
quanto às flutuações à entrada o seu tempo de coe- 
rência tal como é condicionado pela respectiva função 
de auto-correlação. 

A correlação dos valores das sucessivas amostras 
à saída é condicionada imediatamente pela correlação 
dos sucessivos valores das flutuações à entrada na 
medida em que o período de amostragem for compa- 
rável ao tempo de coerência dessas flutuações, e ain- 
da na medida em que o tempo de captação for com- 
parável à constante de tempo do retentor. 

No caso limite de flutuações à entrada de banda 
muito larga, (tempo de coerência desprezável face ao 
período de amostragem), com um retentor ideal em 
que a constante de tempo seja desprezável face ao 
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Fig. 11 


15 kHz 


Função densidade espectral. Resultado experimental obtido para (T,/T)=2,0 e (T,/T)=1,0 


tempo de captação, não há correlação das sucessivas 
amostras e a função densidade espectral à saída é 
do tipo (senx/x)2. A introdução de correlação pelos 
efeitos descritos acima deforma este espectro. 

Verificou-se que os resultados da teoria eram 
confirmados com bastante rigor por medidas feitas 
num dispositivo real, 

Como exemplo duma aplicação do retentor em 
que é fundamental o conhecimento do seu comporta- 
mento em face de ruído, fizeram-se com sucesso ex- 
periências relativas à recuperação da forma de curva 
duma grandeza periódica submersa em ruído, 
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Valores calculados 


Valores medidos 


002 o 02 1 2 2 tr 
Fig. 18 
Densidade espectral na vizinhança da frequência zero 
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Fig. 13 


Montagem para a recuperação da forma de uma curva imersa em ruído 
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Fig. 14 
Circuito RLC usado no ensaio 
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Fig. 15 
Significado do tempo de atraso 7 
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Fig. 16 
Oscilogramas., 
(A) Sinal sem ruído. 
(B) Sinal com ruído. 
(C) Pormenor dum período do sinal (tempo real). 
(D) Forma de curva recuperada a partir do sinal do oscilograma (B), (tempo expandido). 
Comparar com (C) 


5SOmV/div. —— 


iIms/div. —» 


Fig. 1º 
Saída do retentor 
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12. APÊNDICES 
Apêndice - 1. 


Na dedução da expressão da densidade espectral 
duma sucessão de impulsos de amplitudes flutuando 
estocasticamente, seguir-se-á o método da referên- 
cia [5], generalizando-o no que respeita à conside- 
ração da existência de correlação entre as sucessivas 
amplitudes. 

Considere-se para o efeito uma população de 
sucessões finitas de impulsos de duração 'T, cuja for- 
ma é caracterizada pela função g(r) (0<r<T) ea 
amplitude por x. 

A descrição dum membro da população (m) 
como função do tempo será dada por, 


) px ( 

(m mo) ago 

x. (€) = > x, gtt-kT) 
k=—N 


onde N mede a extensão da sucessão que cobrirá um 
intervalo de tempo de extensão (2N+1)T. 

A função de auto-correlação finita deste membro 
da população será então, 


ms [TO dm 
A, () =] A (t) xo (t-=)dt = 


+-N +N die 
(m) qm) [OS - ” 
x x | g(t-p Dg(t-g |-<)dt 
a 


E dao coa Ss 


p=—N q=—N E 


Representando por G(r) a função de auto-corre- 
lação finita do impulso, 


o dos 
G(r) = | g(tigit-:) dt 


e fazendo, 
S=p—q 
tem-se, 
tm) ns ui Qu) (mo) 
Asid= Do D xx GGsT) (AS) 
=—Ns=-p—N 


Passa-se agora ao domínio da frequência, intro- 
duzindo as seguintes transformadas: 


1) Transformada do impulso. 


+D0 
E (um) = ) g(r) exprjwr) dr 
— HO) 


2) Transformada da sua função de autocorre- 
lação 


+20 
f G(r) expí-jor) dr = [(w) [*w) = [U(co)]3 
em) 
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3) Transformada da função de auto-correlação 
da sucessão (m) da população considerada. 


tm) 
N 


+08 (mi) 
Wo ()= A. (7) expí-jor) dr 
-—oo * 


“ 


Esta transformada pode imediatamente escrever- 
-se, partindo de (A-1, 1) e usando a expressão de 2), 


+N  p+N 
Wi =| 
p=-—N s=p—N 


tm) (m 


Mg exp(-jsol) (A-1,2) 


Admitindo que para cada membro da população, 
quando se deixar crescer indefinidamente a sua ex- 
tensão, existe o limite, 


mm) 2 (m), 
E O cerne NY ) 
(0) Ein (QN41) Tn (03) 


P 
ele traduzirá a densidade espectral correspondente 
a esse membro, expressa por Hz de largura de banda, 
isto é de tal modo que o valor quadrático médio da 
sucessão x(”) (t) seja dado por, 


(tm)? "20 (uu) do 
(xo ) já |, ua 


Se se admitir a estacionaridade da estatística das 


sucessões de amplitudes x * o somatório interno 


da expressão (A-1,2) será independente de p, sendo 
o resultado do somatório externo em p simplesmente 
o resultado do produto do somatório interno em s 
pelo número de parcelas em p, isto é, (2N+1), pelo 
que se obtém no limite N -—s> 0. 


2 es 
tm) ,. tm) (m) ma 
to  ESOAES 7) 3 - 
P “(u) T E (eo)| + x. Xs exp(-jsol) 
S = O 


Finalmente, fazendo a média sobre a população 
e representando por X o valor expectável do pro- 


duto XX -, que só dependerá de s, 


x. = CX X > 


tem-se para valor expectável da densidade espectral, 


+00 
> ha exp(-jsol) 


$S =——CO 


<P(o)>= = EF (m)|? 


ou ainda, dado que X será par na sua dependência 
de s, 
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+ 


9 OS 
<P(u)>= | F(o)|t | Xo +? 2 X cos(suT) 
s= 


Admitindo que se trata de fenómenos estocásti. 
cos ergódicos, as médias sobre a população podem 
ser igualmente encaradas como médias ao longo do 
tempo relativas a um qualquer membro da população. 


Apêndice - 2. 


Considere-se um filtro passa baixo de função de 
transferência racional G(s), excitado por flutuações 
estocásticas de espectro branco caracterizado pela 
densidade espectral constante, independente de «w, nº 
(em V2Hz -1). 

Nestas circunstâncias a função de auto-correla- 
ção das flutuações à saída será dada por, 


42 e ==s 


joo ds 
3 o G(s)G(-s) exp(sr) e 


A função F(r) é par e pode ser obtida no domi- 
nio 7>0 completando o eixo imaginário, ao longo do 
qual está definido o integral, com a semi-circunfe- 
rência no infinito do semi-plano Re (s)<0, que nestas 
circunstâncias não contribui para o integral esten- 
dido ao contorno fechado assim formado. 

O resultado pode assim ser expresso imediata- 
mente em termos dos resíduos da função integranda 
nos polos abarcados pelo contorno, 

No caso presente, trata-se de dois polos reais 
distintos correspondentes às constantes de tempo T, 
eT.. 


l 1 


E 3 UM ==. E — 
(s)G( s) (1+sTo) (15 sEyJ( - so) (1—sT,) 
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Calculados os resíduos nos polos do primeiro fac- 
tor obtém-se, 


2 1 1 
F(y)=—| —— ——*— exp(— |el/To) + 

A) 7 Tº aê q” 

t 
y (A-2, 1) 
1 Ea 
—— -— |xl/T 
1 0 


Se nesta expressão fizermos T, — 0, obtém-se a 
função de auto-correlação das flutuações y(t) à en- 
trada do retentor, contando pois com o efeito apenas 
do filtro que precede o retentor, e que será então, 


n2 na 
Y(7) Es 4T exp(—!r|/ l 1) 
1 


O valor quadrático médio das flutuações à en- 
trada do retentor obtém-se tomando esta função 
para 7=0, 


2 q? 
Tm UU) = — 
ef 4T, 


o que permite substituir em (A-2,1) nº pelo seu valor 


2 
expresso em termos de Y ' 
“ 


r2 2 
A 21 


obtendo-se assim o resultado expresso por (6.2) e 
(6.3) no texto. 
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Teorias de cascas envolvendo momentos 
“ * 
de ordem superior à primeira É) 
EDUARDO ROMANO DE ARANTES E OLIVEIRA (') 
RESUMO SUMMARY 


Tendo o desenvolvimento extraordinário dos mé- 
todos numéricos, e em especial do método dos ele- 
mentos finitos, despertado interesse pela construção 
de novos modelos estruturais, é objectivo do trabalho 
apresentar e comparar duas teorias de cascas nas 
quais as tensões generalizadas são as resultantes e os 
momentos de qualquer ordem das tensões distribuidas 
na espessura. Na base do trabalho está a considera- 
ção de que o conhecimento dos momentos de ordem 
superior à primeira permite mais precisa determina- 
ção da distribuição das tensões que a que pode ser 
obtida de uma teoria de tipo convencional. 


1 — O INTERESSE DE TEORIAS ENVOLVENDO MOMEN- 
TOS DE ORDEM SUPERIOR A PRIMEIRA 


Uma função que se distribua uniformemente num 
intervalo pode ser completamente definida nesse in- 
tervalo pelo seu valor médio, 

O valor médio (ou momento de ordem zero) e o 
momento de primeira ordem relativamente a um dado 
ponto chegam por outro lado para definir uma função 
que se distribua linearmente. 

De uma maneira geral, uma distribuição polino- 
mial de grau N, como qualquer distribuição depen- 
dente de N+4-1 parâmetros, é completamente definida 
pelo conjunto dos seus momentos de ordem não menor 
que zero e não maior do que N. E um conhecido teo- 
rema afirma que qualquer distribuição limitada é 
univocamente determinada pela sucessão de todos os 
seus momentos. [1] 

Ora, sendo as tensões generalizadas, nas teorias 
convencionais das peças laminares (ou cascas), re- 
sultantes e momentos de primeira ordem das tensões 
efectivas (”), as referidas teorias só poderão permitir 


Considering the interest of new structural models 
which the expansion of the numerical techniques of 
analysis, like the finite element method, very much 
increased, the aim of the paper is the presenta- 
tion and comparison of two theories of shells in which 
the generalized stresses are the resultants and higher 
order moments of the stresses distributed in the 
thickness. The basic idea of the paper is the fact that 
the knowledge of higher order moments permits a 
more precise determination of the stress distribution 
that the one which can be provided by a theory of 
a more convencional kind. 


o cálculo destas tensões se se lhes associar uma hi- 
pótese que faça depender de dois únicos parâmetros 
a sua distribuição na espessura, 

Tal é por exemplo a que lhes admite uma distri- 
buição linear, hipótese esta que se aproxima excelen- 
temente da realidade no caso das cascas finas, isto é, 
daquelas em que a espessura é muito pequena em 
comparação com os raios de curvatura ou com outras 
dimensões características em jogo, mas já não no 
caso das cascas espessas, para a análise das quais 
várias teorias [2] têm sido apresentadas que podem 
ser olhadas como aperfeiçoamentos da teoria das 
cascas finas na forma que Love lhe deu. 

Tais complicações não devem porém fazer es- 
quecer o ponto fundamental de que qualquer hipótese 
que faça depender de dois únicos parâmetros a dis- 
tribuição das tensões na espessura não pode repre- 
sentar mais que uma simulação muito grosseira da 
realidade. 

Por outras palavras, só teorias envolvendo mo- 
mentos de ordem superior à primeira poderão então 
conduzir a resultados satisfatórios, isto é, a resulta- 


(*) Comunicação apresentada ao 1.º Congresso Nacional de Mecânica Teórica e Aplicada. 

(1) Professor Catedrático de Engenharia Civil, investigador do CEMUL, 

(*) Exceptuam-se as tensões tangenciais transversais cuja distribuição (parabólica anulando-se nas extremidades do 
intervalo) depende de um único parâmetro, e às quais se associam portanto somente as resultantes, vulgarmente designa- 


das por esforços transversos. 
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dos suficientemente aproximados dos fornecidos pela 
teoria tridimensional, a qual não fixa em princípio, 
como é sabido, qualquer limitação à variação das 
grandezas na direcção perpendicular à superfície 
média. 

O presente trabalho tem por objectivo apresentar 
e comparar duas dessas teorias com base em relató- 
rios que, embora elaborados pelo autor há já alguns 
anos, nunca tiveram mais que uma difusão res- 
trita (3, 4]. 


Observa-se que o que vai dizer-se sobre cascas 
poderia ser dito sobre barras. A adaptação não seria 
difícil, e só para não complicar mais o trabalho nada 
se dirá sobre o assunto. Não se deixa de mencionar, 
no entanto, que uma das aplicações mais interessan- 
tes de uma teoria com momentos de ordem superior 
à primeira parece poder fazer-se na análise das cha- 
madas vigas-paredes [5], Não houve infelizmente tem- 
po para a discutir na presente comunicação, mas o 
autor não deixará de a incluir nas suas preocupações 
futuras. 


2 — SISTEMA DE COORDENADAS 


Serão usadas coordenadas «,:, E tais que a 
superfície coordenada £ =0 coincida com a superfí- 
cie média S e qualquer linha correspondente a valores 
constantes de £, e É, seja recta e perpendicular a S. 


Simbolizam-se com h, (E, E,» E) 08 factores de 
escala (parâmetros de Lamé) associados com as 
coordenadas é, e faz-se h, constante e igual à uni- 


dade. 


Supõe-se ainda que as linhas coordenadas per- 
tencentes à superfície média são linhas de curvatura 
principal dessa superfície, o que tem por consequência 
que o sistema de coordenadas seja ortogonal em 
todo o espaço e os versores tangentes às linhas coor- 
denadas deixem de depender da coordenada E, isto é, 


= O (1=1,2,3) (2.1) 


a 
RAR 


Simbolizando com ho, ho, Ro e Ro os factores 


de escala e raios de curvatura principais em pontos 
da superfície média, tem-se 


x (a = 1, 2) (2.2) 


Supõe-se que os índices gregos podem tomar os 
valores 1 e 2, podendo os latinos tomar os valores 
1,263. 

Observa-se que as grandezas referidas à super- 
fície média (assinaladas com o índice *) poderão ser 
consideradas, já que correspondem a um valor fixo 
de é, (£,=0), como funções exclusivas de £ e £. 
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Represente-se agora por À O cociente entre o 
factor de escala ho (E, £,» $,) num ponto arbitrário 
da casca e o factor de escala h, (E,,E,) no ponto 
da superfície média com as mesmas coordenadas E 
Demonstra-se[6] que 


Ea 


Rº 


EA 


(2.3) 


à mta 
a 


A fronteira da casca, considerada como um cor- 
po tridimensional é composta por duas partes. À pri- 
meira, =, é formada por porções de duas superfícies 
com equações E =-+12e(;,,5%) em que e repre- 
senta a espessura. A segunda, E», é-o por uma ou mais 
superfícies com equações do tipo f (£,, £,)=0. 

Designe-se por I a intersecção de s, com 8. 

Qualquer ponto de XE, pode ser referido a um 
sistema de coordenadas £, £., sendo as linhas coor- 
denadas correspondentes a £, rectas e normais a T. 

Representa-se por à a relação entre os valores 
do factor de escala h associado com a coordenada £ 
num determinado ponto de =, e o valor do mesmo 
factor de escala no ponto de T com a mesma coor- 
denada é. Tem-se pois 


h (5,54) 
hº (E) (2.4) 
à desempenha evidentemente o mesmo papel em E, 
que as grandezas A nas superfícies coordenadas 


É Pa constante. 


3 — UMA PRIMEIRA TEORIA 


Qualquer modelo estrutural deve apresentar o 
mesmo esquema que a teoria tridimensional, [7] isto 
é, três grupos de equações (equilíbrio, deformações- 
-deslocamentos, tensões-deformações) envolvendo três 
espécies de grandezas (tensões-deformações, desloca- 
mentos) que se denominarão generalizadas (no sen- 
tido de Lagrange) por oposição às efectivas que apa- 
recem no modelo tridimensional. Supõe-se por outro 
lado que se verifica o princípio do trabalho, o qual, 
combinado com as equações de equilíbrio, fornece as 
equações deformações-deslocamentos, e, com as equa- 
ções deformações-deslocamentos, fornece as equações 
de equilíbrio. 

A melhor maneira de estabelecer uma teoria en- 
volvendo momentos de ordem superior à primeira 
será justamente começar por estabelecer as equações 
de equilíbrio e deduzir seguidamente as relações de- 
formações-deslocamentos recorrendo ao princípio do 
trabalho. [6] 

O primeiro problema que se põe é pois o de esta- 
belecer as equações de equilíbrio, as quais deverão 
necessariamente ser exactas ou, mais precisamente, 
deduzidas das equações de equilíbrio da teoria tri- 
“dimensional. 
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Estas, em coordenadas curvilíneas ortogonais, 
são 


(h, h; 3)n+ (ha h, co + (h, ho 0a)»a + h, h, h,X — O 


A utahoataAa+ Y=O (3.2) 
em que os vectores o, são os vectores das tensões 
efectivas actuantes nas faces normais aos versores a, 
e os vectores X e Y representam as densidades das 


forças e momentos exteriores distribuídos no volume 
do corpo. 
No presente capítulo supor-se-á Y=0, o que equi- 


vale a admitir que o tensor das tensões é simétrico, 
isto é, que c,,=0,;. 

Atendendo ao que ficou escrito em 2, a equação 
(3.1) pode transformar-se em 


(h; ha Sim + (hº À, Saba (hy h; à da Saya + 


+h h,44X=O (3.3) 


Multipliquem-se agora externamente ambos os 


membros de (3.3) por E, a,e integrem-se seguida- 


mente relativamente a &, entre sy 047 E, Obtém-se 


e), 


n 
f és a À b (hz “ A + (hyhoa, à Saba + 


-€/, 


* By ça aXJau—o (3.4) 


em que q é igual a 2 quando q for igual a 1 e igual 


a 1 quando « for igual a 2. 
Transformando [3] a equação (3.4) e supondo 
ns 1, obtém-se 


Di Lina ANO AMP bia, ANDO )+ 
3 


tn) 


+h/hím-0O (3.5) 
em que 
(m) (mn) 
NGC (3.6) 
a Va 
sendo 
e/2 
(mn) | n 
N = A oc dê, (3.7) 
ab Z ap 
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e/2 


(0) [a 
N, =| E (l-+na) Az Tag ds (3.8) 


[-4 
-e/2 


O vector mn), função das forças exteriores 


actuantes na casca, tem a expressão 


e/2 
(1) n a 
m =a, A | [enusa 
-e/s 
eis 
n e 2 e 2 
És du Qu:6 | = —. | 
+fó REV + (E) 


«8/2 


(3.9) 


a qual, evidentemente, se simplifica se a espessura 
for constante. Observa-se que m(n) não tem compo- 


nente paralela a a ,. 


A equação (3.5) é válida, como se disse, para 
ns 1. Para n=0, pode estabelecer-se facilmente a 


equação 


O (0) 0,0 fo) 
2 pn Ja A h, p=O 


(3.10) 


na qual 
e/2 
(0) 
PO = A NX dis + 
-e/2 
o/2 
Cn 1º /€u 1? 

Es n2/1+ (E)A(G2) (3,11) 

-p/2 


Observa-se que tanto a equação (3.10) como a 
equação (3.5), particularizada para n=1, aparecem 
na teoria convencional. 

Estabelecidas as equações, podemos agora indicar 
as tensões generalizadas que intervêm na teoria. São 


elas as resultantes e momentos JN vo? que aparecem 


nas equações de equilíbrio. 
Sendo N o máximo valor de n, as equações de 
equilíbrio envolvem todas as grandezas N pr desde 


n=0 até n=N e todas as grandezas N (1) desde 
n=0 até n=N-1. As grandezas N (x) não intervêm 


nas equações e não serão pois tomadas para tensões 
generalizadas. 

As tensões generalizadas são pois 6 N+4 e as 
equações de equilíbrio 3 (N+1). 
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Na teoria convencional, tem-se N=1, pelo que 
as tensões são 10 e 6 as equações de equilíbrio. 

As equações de equilíbrio a satisfazer no con- 
torno TF, ou sobre qualquer linha traçada sobre S, 
podem ser obtidas procedendo com a equação 


(3.12) 


do mesmo modo que com a equação (3.3). Obtém-se 


(uu) (nm) 


do A adm E (3.13) 
para ns1e 
para n=0, sendo 
o/2 
(mn) 
: a, À f Ago dE, (3.15) 
-9/2 
e/2 
(0) |º 
N = | à a dE, (3.16) 
“ei? 


As equações deformações-deslocamentos podem 
agora ser estabelecidas com toda a facilidade [3] uti- 
lizando o princípio do trabalho cuja expressão é 


N f 
im) (n) (0) (0) 

pô 2 ENA as | Pp U + 

n=0 a « ' 
s s 
N 
(nm) Am) (0) — to) 

+ m.v as + | N.U + 

n=1 á di: 

r 


(3.17) 


As tensões N'v supõem-se equilibradas pelos 
campos de dorgua: extrrioen definidos pelos vectores 
p (o) é m tn) em S e Nt é M tr em T. 

As deformações E tu definem-se como os multi- 
plicadores das tensões NU vu) na expressão do traba- 


lho das forças internas, e os deslocamentos como 
multiplicadores das forças na do trabalho das forças 
externas. 
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Introduzindo as equações de equilíbrio em (3.17) 
e aplicando o teorema da divergência, obtém-se 


N + 
(n) (tn) (0) 
2 2 Rê E, dS-D [nt 
a n=0, x ; 
s s 
to) N= 
U (m) 
ms as AN a dS + » N, 2 
hz n=1 
s 
(u) 
ta tm «+ 1) 
Aa — a, |dS+ 
o a é es 
ha 


(3.18) 


e, como esta equação deve verificar-se para valores 
contínuos mas arbitrários das tensões N (n) , obtêm- 


-se as seguintes 6N +4 equações 


U, 
ed ua o ve ” (3.19) 
o (0) 32 a 
ha 
vi 
(m) 
Ea ta "a Aid (3.20) 
e (0) - 3 - e X 
x 
(N) 
(3) Via 
Esp — TN As Ag (3.21) 
h 


que relacionam 6N+4 deformações generalizadas 
com 3(N+1) deslocamentos generalizados, Para N=1, 
obtém-se as conhecidas equações da teoria conven- 
cional. 

Quanto às equações tensões-deformações, elas 
têm a forma 


(3.22) 


em que a função wº (N pr ), densidade bidimen- 


sional da energia complementar, pode ser determi- 
nada a partir de 
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ei 
MO) = | w E o | xrade (3.98) 
2 mu : s lj a a 
" -e/2 


em que W, é a densidade tridimensional. 

As equações tensões-deformações (3.22) poderão 
ser desenvolvidas somente se for introduzida uma 
hipótese que permita exprimir as tensões efectivas 
em termos das generalizadas. 

Finalmente, as condições de fronteira são 3+2N, 
a saber 


NÍO vo = N69) (3.24) 
(n) o (nm) 
a AN vVv=M (3.25) 
sobre um troço T, deT e 
(0) to) 
U=U (3.26) 
( cn) 
Rs (3.97) 
a “a 


sobre o troço complementar, 1, . 


4 — UMA SEGUNDA TEORIA 


A teoria bidimensional que acaba de ser apre- 
sentada representa aparentemente uma aceitável ge- 
neralização da teoria convencional (correspondente 
a N=1). 

Uma crítica cuidadosa não pode deixar no en- 
tanto de lhe encontrar o grave defeito de a sucessão 
das soluções geradas por sucessivas teorias bidimen- 
sionais correspondentes a valores crescentes de N 
não tender para a solução correspondente à teoria 
tridimensional quando N tende para infinito. 

Para concluir que assim é basta notar que, dado 
que a componente X, da densidade das forças de 
massa não intervém nos vectores m("), a equação 


(3.5) não envolve esta componente. Por este motivo, 
as soluções bidimensionais não dependem do modo 
como X, se distribui na espessura. Ora, como tal não 
acontece com a solução tridimensional, segue-se que 
não pode esperar-se atingir a solução exacta fazendo 
tender N para infinito. 

Igual conclusão se poderia tirar do facto de a 
equação (3.13) não envolver a componente o, do 
vector o. 


Entre as tensões generalizadas não aparecem por 
outro lado momentos das tensões o,., o que torna 
impossível o cálculo destas em termos daquelas. 
Observa-se que isto já acontece na teoria convencio- 
nal que também não fornece qualquer informação 
sobre as tensões o,.. 
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Na teoria que seguidamente se apresenta elimi- 
nam-se as imperfeições apontadas na primeira. 

Parte-se como anteriormente das equações de 
equilíbrio (3.3). 

Multiplicando ambos os membros por £,* e in- 
tegrando obtém-se: 


e/2 
[DCE + 


x 


lg * 
e/2 
di CIÉRT A 8 di + 
+ 4 E EO 2 3 dos 5 
-e/2 
e/2 
o n 
+ hh lxxét Xdé = O (4.1) 
ê «8 O E Mus 5 
-e/2 


equação esta que, transformada, dá 

o (un) o o n-1) (n) 
> (h N )-h ME q TES p —O (4.9) 
p ã E ai 1 2 Ea 1 2 - 


desde que n 51. 
Para n=0, tem-se, em vez de (4.2), 


po o) (0) 
5 N) +hhp=0 (4.3) 
a aaa LB E 
sendo 
e/2 
(n) n 
Pp = Rx t hã dt + 
O E Má 3 
-e/2 
a en e» 42 |º12 
+[ Eras fo +(3) | (4.4) 
3,3 & h, / h, -e/2 
por qualquer valor de n 5 0. 
As componentes dos vectores Na NO não 
“ meé 


devem ser consideradas como completamente inde- 
pendentes. Com efeito, supondo à Cet a 0, Mimi AR 
e, portanto, obtém-se 


(n-1) 1) NaS 
E 


o) 4. 
a o (4.6) 
Por outro lado, como 6 = T, 
Ntº) Ntt) 
(n-1) (n-1) 12 91 (n=1) 
Ee DE aj e fr E e E | O (A) 
12 21 R 1 R a 5 
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Dn 
E | “it SAL (4.7) 
«0/2 


Introduzindo (4.5) em (4.2), obtém-se 


(0) (0) 0.0 (n-1) 
> (1 N ).. -h hn ae, + 


x 
J 


o o (n) 
3 | 2 ia 


As equações (4.3), (4.6) e 4.7), consideradas para 
valores de n não superiores a N, formam o conjunto 
das equações de equilíbrio da teoria. O seu número 
total é 4N-+3. 

As tensões generalizadas são todas as resultan- 
tes e momentos envolvidos nas equações de equilí- 
brio. Como estas contêm todos os vectores N (?) 


para 04n<N, tem-se 6 (N+1) tensões generaliza- 


(n) (n) 
das Nai. Aparecem por outro lado N tensões N., 
para 0<n<N-1, o que perfaz um total de 6 (N+1) + 
+N=7N+6. 

As equações de equilíbrio a satisfazer sobre T 

ou qualquer outra linha traçada sobre S, 
(nm) n) 

N 


nal 


x 


É au N (4.9) 


obtêm-se a partir de (3.1) multiplicando por À £,º e 


integrando, Observa-se que 


eiz 


(m) n 
N - [41 a dé 
E a 3 


-e/2 


(4.10) 


A utilização do princípio dos trabalhos virtuais[4] 
conduz às seguintes 7N+6 equações deformações- 
“deslocamentos 


(0) - (1) 0) 
= "+U a-V aaa (4.11) 
- & E a - q 3 “a a 
(n) 
(n) e x (n+4+1) 
x ha a E 
[4 
(n-1) 
(0) Vs, 
—a A a (v + — ) para0<n<N (4.12) 
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(x) (N) 
AS) a U, di 
E = — + N = & — 3 A q 2 (4.123) 


(n) 


(n+1) 


nas quais V.(") representa a componente de desloca- 
mento associada com m, (º). 
As 3 (N+1) componentes dos vectores U(") e as 


N grandezas V.("), multiplicadores dos parâmetros 
definidores das forças exteriores na expressão do 
trabalho, são os deslocamentos generalizados. 

As 7N+6 deformações generalizadas, multipli- 
cadores das tensões, são por sua vez componentes 
dos vectores Un) e as grandezas V,tn). 

As equações tensões-deformações são as equações 
(3.22) e ainda 


po E Wo, (4.15) 
E (| 
“iNaa 


Quanto às condições de fronteira, são 3 (N+1) 
em cada ponto de T, a saber, 


(m) —(n) 
vu =N em TF, (4.16) 
a & 
(n) —(n) 
U = N em T, (4.17) 


5 — CONCLUSÕES 


Foram apresentadas duas teorias. 

Em ambas, N representa a ordem máxima dos 
momentos intervenientes. 

A primeira contém 3 (N+1) equações de equi- 
líbrio, 6 N+4 equações tensões-deformações e 6N+4 
equações deformações-deslocamentos, envolvendo 
6 N+4 deformações e 3 (N+1) deslocamentos. Ao 
todo, 15 N+11 equações e 15 N+411 incógnitas. 

A teoria convencional aparece como caso parti- 
cular desta teoria, correspondente a N=1. As solu- 
ções fornecidas pela teoria não convergem porém para 
a solução tridimensional quando N tende para infi- 
nito. 

A segunda teoria contém 4N+3 equações de 
equilíbrio, 7N-+-6 equações tensões-deformações e 
7N+6 equações deformações-deslocamentos, envol- 
vendo 7 N+6 tensões, 7N+6 deformações e 4N+3 
deslocamentos. Ao todo, 18 N +15 equações e incógni- 
tas. 
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Esta teoria já não fornece a teoria convencional, 
uma vez particularizada para N=1, Envolve por 
outro lado mais incógnitas e equações que a primeira. 

Estes inconvenientes são compensados pelas van- 
tagens de a teoria dar informação sobre a distribui- 
ção de o,, na espessura e, sobretudo, de as sucessões 
de soluções correspondentes a valores indefinidamen- 
te crescentes de N convergirem para as soluções 
geradas pela teoria tridimensional. 

Como os valores de N devem manter-se pequenos 
nas eventuais aplicações, é possível que esta defi- 
ciência de convergência não tenha consequências prá- 
ticas apreciáveis, Só a aplicação à resolução de pro- 
blemas concretos permitirá tirar conclusões, 
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Aplicação da solução do problema de Mindlin na resolução 
do problema da estaca solicitada por uma força horizontal 


SUMÁRIO 


Com base na solução do problema de Mindlin, 
que considera o campo elástico gerado por uma força 
concentrada actuante num ponto interior dum semi- 
-espaço, do qual se deduzem por integração os deslo- 
camentos devidos a uma carga distribuida linear- 
mente, apresenta-se um método para a resolução do 
problema da estaca cravada num solo homogéneo 
solicitada por uma força horizontal. 


1 — INTRODUÇÃO 


A resolução do problema da estaca solicitada por 
uma força horizontal no topo tem sido proposta por 
vários métodos e em especial pelo dos elementos fini- 
tos. Neste, o número de elementos a empregar, para 
uma razoável aproximação de resultados, é geral- 
mente elevado, O método que se apresenta nesta co- 
municação, mais exacto embora menos geral, permite 
a comprovação, no caso de meios homogéneos, dos 
resultados obtidos pelo método dos elementos finitos. 


Considerando-se válidas as hipóteses da lineari- 
dade física e geométrica, introduziu-se a hipótese 
simplificativa de a distribuição contínua correspon- 
dente às forças de contacto entre a estaca e o solo 
ser horizontal e paralela à força actuante no topo 
da estaca. A distribuição referida foi discretizada 
num conjunto de cargas triangulares elementares, im- 
pondo-se então a compatibilização, segundo a direc- 
ção da força, dos deslocamentos da estaca, conside- 
rada como peça linear, e do terreno, num número 
finito de pontos do eixo da estaca. Aplicou-se por- 
tanto o método dos esforços. 


A determinação dos deslocamentos produzidos pe- 
las cargas triangulares elementares implicou a reso- 
lução do problema da determinação dos efeitos pro- 
duzidos por uma carga distribuida linearmente numa 
zona interior dum semi-espaço elástico. 


ANTÔNIO JOSÉ REIS (') 
EDUARDO BORGES PIRES (*) 


SUMMARY 


The solution of Mindlin's problem, which consi- 
ders the elastico field generated by a force concen- 
trated in a point located within an elastic half-space, 
from which the effects of forces distributed on lines 
can be generated, is used for the solution of the 
problem of the pile acted by a horizontal force. 


2 — DESLOCAMENTOS PRODUZIDOS POR UMA CAR- 
GA LINEARMENTE DISTRIBUIDA NO INTERIOR 
DUM SEMI-ESPAÇO ELÁSTICO 


21 — O problema de Mindlin 


Designemos por z o eixo normal à superfície do 
semi-espaço que contém o ponto de aplicação da for- 
ça, a qual actua paralelamente à fronteira, segundo 
o eixo dos xx, fig. 1. O vector deslocamento num 
ponto genérico do semi-espaço, u=ue,+v e,twe, 
exprime-se [1] por: 


1 
———— v (r.B-+8) (2.1.1) 


41 —»v) 
em que » é o coeficiente de Poísson do meio, Be 8 


as funções de Papkovitch, y o operador vectorial 
g ê q 


V=—"— € Re Dm 
há ext E Sy E duda 


erovector r=Xxe +yec +ze 
al EK “Bs 
Desprezando as forças de massa as funções B e A 


verificam a equação de Laplace ou seja: 


vB =0 e v2B = O (2.1.2) 


(*) Comunicação apresentada ao 1.º Congresso Nacional de Mecânica Teórica e Aplicada. 
() Eng.º Civil, Assistente do I. 8. T., investigador do CEMUL. 
(*) Eng.º Civil, Assistente do I. S. T., investigador do CEMUL. 
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Com base nas equações (2.1.2), Mindlin [1] de- 
terminou directamente, por aplicação da teoria do 
potencial, as funções de Papkovitch para o problema 
em questão: 


(2.1.3) 


x(1—2v)? 


P[ x1-2,) ex. 
4 , e E, 
2:G | R(R+2+0) R 
| (2.1.4) 
Ratz+c 


SE 
Pe “Tal RR ++) 


em que G é o módulo de distorção do meio, c a coor- 
denada z do ponto de aplicação da carga e R, e R, 
as funções: 


R=[xº+y?!+(2z-— cc) 
(2.1.5) 


R = [2+y'+(z+ep]” 


Os deslocamentos (2.1.1) são portanto infinitos 
no ponto de aplicação da carga. 


2.2-— O problema da carga distribuida linearmente 
numa zona interior dum semi-espaço elástico. 


Atendendo às hipóteses simplificativas introduzi- 
das em 1, determinaram-se apenas as componentes 
u do vector deslocamento u=ue,tve,twe ; Sobre 


o eixo do zz. Deste modo designando por u, a com- 
ponente u do deslocamento devido à carga concen- 
trada P e atendendo à equação (2.1.1) temos: 


: :B, :B, 
4(1— v)| = ix ix + 


(2.2.1) 


Atendendo às equações (2.2.1), (2.1.3-5) e à fig. 1 
os deslocamentos u, em x=0 serão: 


A ) - a (m ES 
! O di R, + 2) — 
P VA 1-2 y E. 2 
4( 1 — v ) 27G RAR, + Zz + Cc) Rº mas 
2 
1—2y 1-2y c 


(2.2.9) 


— WG | Rfzte * RiR,4240) 
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Os deslocamentos u em x=y=0 devidos a uma 
carga linearmente distribuída p=p, +p,c em (— a < 
<y<a; ec <z<c,) serão dados por: 


n 
(Dececo = AG +Sy) — 7  feBra a og 
-SJ4DI(I-?)S+Si) (228) 
em que: 
chi | 124 
-— O | B=5E | D=-E 
= l6rG(1I— py) 27G 27G 
(2.2.4) 
—DT Pd tas) 
o [s2+ (z-e)']t! 
“e 
— MP qd (220) 
ul Cs [st + (z+c)2)!!? 
S, = IN ii 
dad q (Ste (z+e)]to ([ Hz + e) ty z+e) 
de ds 
(2.2.7) 
—PotPO) 4 q, (228) 
o [83 + (z+erPp 
E. —— TA dcds (2.2.9) 


Hz +ejjtP+z+e 


= ç Es 
Tri Pe) 


dc ds 


(2.2.10) 


No cálculo dos integrais (2.2.5-10) fomos condu- 
zidos frequentemente a integrais impróprios na variá- 
vel s, pois as funções integrandas apresentam descon- 
tinuidades infinitas para s=0. O cálculo processou-se 
então geralmente através de limites o que, atendendo 
ao facto das funções integrandas serem funções pa- 
res, relativamente à variável s, se efectuou de acordo 
com a expressão geral: 


“a “es 
NE sã lim | J é (se) de ds, i= 1,2,.. 55 
E-» () E 
(2.2.11) 


em que q, (sc) representam as funções integradas 
em (2.2.5-10). 

Apresentam-se seguidamente os resultados dos 
cálculos efectuados, tendo-se designado por aj=|2-C,|, 
Bi=z+ce, e y,=c-z, com i=1,2. 
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Ps ipaKi 
P cao campo 
8533 Kgterê 


go »10em 
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220 


240 
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FO 8 


utmem) 


— ulz)- corgo fneormente distnbuida 
-—- WAZ)- carga com. estoticamente equivalente 
esc dos uu - 112 


esc dos 11-3H% 
esc de fo<as- tem O 100K1 


SECÇÃO DA ESTACA 


«bs Da 


(Do * Pyc) dc ds 


dE 4 


pn bem 


tu 


FO 1 
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z1<c, >S, = 2a E (e E) (Po + P4Z) + 
et (ate) 


1/2 
placa 
tá E (et) (2Po * Ph) — 
%, 
, ' gua E 1/2 
— à, | log (Meta) (2po + P;) + 


fia [e (use)! | nam 


[= a ta) 
z=e>S,=º?a(p+polog | -—+ 20º 1+ 
a 


a+(a?+ a) e 
+aiPo + Pájllog| —+—— 2 14 
ã 


a aa] 


B<Zz<eo— S, = 28x 


(20) [es (0) 


(2.2.13) 


>< log 


é 1/2 
+( a? fa 
>x<(Po + Py2) + as(2po + Pisa) X log Ema 
&a 


o Wi 
+ «(2po + p,8,) log [+ | + 
+ap, (nº + a) jr (a: E, "| (2.2.14) 


+ (80) 


(+) 


al (a? 4a ” 
aj(2Po + PiSg)log [eta + «(2po+P,8,)>x< 
Z, 


eos | HLnt o) +a)" Jemmtto ss - 
(ee 


(2.2.15) 


a, + (a: + a)o 
z<€ — 5, = la o (Po + PiZ) + 


a+(u+a |, 


+ a(2Po + P;2,) os? 


+ ap, [a —(aº + a) | (2.2.16) 
21682112 

S,=2a(Po - Zp,) log di +B(2Po+yaP1)>< 
8, +(a? 185) 1/2 


2 


aa E ter 
(1246) 


218 
uq (a: + B ) 
>< log sea —B,(2Po + Px 


(2,2.17) 


ê a + (a + 6) 


S, = 2(Po o P,Z) log qn ee 
3 É a + (a +86)! 


Aids Gi E (aº + 8) 


Life sy? 
B + (2º + s)" FR +8) ii 


>x(2po + YaP;) — (a + By” (2Po + YP) — 


a, [2po + pi(es + e] (2.218) 


2 41/ 
218 )º 
S =2 jor-tno log É | + 


1 a-+ (aº = 6) 


214/ 
B+(a +86): 
-+ap, log [Eta] + 


p+(2 +86)! 


+ (Po pa] (2 +8)1— 


CR] 


(2.2.19) 
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Ba + (at “je 8 ja 
S;=a(po—P;z) log] —— 20 + 


+ 8a(2p, + YoP;) log | 


i+ (a 4 8) 
— B(2Po + nP;) log | — + 1 + 


+ (po— pa) [8 88, (a +68) 4 
enforep)eRafleso)P (e m)- 


(808) (28) a(o! = 8)] aaa 


Se = — 2z(Po — Pp;Z) log 


pr (atg)! 
+ (Po — 2zp;) log | — 13 


9 1/2 1/2 
—E (po — apo >| (2 +85) 0 (at 8) 0 

a 2, 
Va 


f 


cd) (a + E)! o 8 — s)| + Pi s(u o)! - 


a (ur + “go 8 - 8º) aa (uº E 5)" 


en 


A verificação das expressões foi efectuada tendo 
em consideração que estas devem tender para as 
correspondentes expressões contidas em (2.2.2.), 
quando a carga distribuída p=p, +p,c tende para 
uma carga concentrada P aplicada em z=c, ou seja: 


a fáio e) [+ (Po+ ep) [o(a +8 


(2.2.21) 


, Pp 
= 8 


com 1=1, 2,...,6 (2.2.22) 


em que S. são as expressões correspondentes à 
carga concentrada. 
Além da verificação (2.2.22) efectuou-se ainda 


uma verificação numérica. 
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Efectivamente os deslocamentos u em x=y=0 
correspondentes a uma carga distribuida numa pe- 
quena zona interior dum semi-espaço devem ter uma 
distribuição semelhante aos produzidos pela resul- 
tante estaticamente equivalente, a menos duma vizi- 
nhança da zona de aplicação da carga. Tal é asse- 
gurado pelo princípio de Saint Venant. Os resultados 
desta verificação apresentam-se na fig. 2, 

O deslocamento (u)x=y=0 fica caracterizado 
num dado meio (y, G) pelos seguintes parâmetros 
Cc,» €C Pp P, €Zz0u seja: 


(u) =uU(Cc, C PyP,p2Z) (2.2.23) 


z=y=0 


3 — FORMULAÇÃO DO PROBLEMA ESTRUTURAL PELO 
MÉTODO DOS ESFORÇOS 


3.1 — Generalidades. Vector dos deslocamentos no 
sistema base 


A acção do terreno sobre a estaca foi simulada 
por sobreposição de cargas triangulares elementares 
conforme dissemos em 1, Seja então o sistema base 
representado na fig. 3 b) no qual a carga F é equili- 
brada pelas cargas triangulares representadas, de 
densidades máximas respectivamente: 


F(2l + z,) Fz, 
PaHiO a(i=z, ) (Ql+z, 2) 
(3.1.1) 


Po - 


“za(a-+z,- z,) 


O deslocamento u , no nó i do sistema base, 
decompõe-se em duas parcelas. A primeira u",, de- 
vido a um deslocamento de corpo rígido da estaca 
o qual é originado pelas deformações dos «apoios 
fictícios» extremos, respectivamente Voo e u”, (ada) 
A segunda devido à deformação de flexão produ- 
zida pelas reacções. Temos então, 


= Us mé (3.1.2) 
com 
quim U 
uU =u du e ,, (3.1.8) 
oi 00 l 


Os deslocamentos Vo são os desloca- 


E Us(n+1) 
mentos produzidos no terreno pela acção simultânea 
de cargas iguais mas de sentido oposto às reacções 
de apoio, os quais são obtidos pela equação (2.2.3). 

Designemos por U o vector dos deslocamentos 


no sistema base, 4 = (U,,) 1=1, 2, ..., n. 


125 


126 


so 4 


um 

TE 
1/0 INRANNARE 
E Naa 


Nu 


Es + 


Jj 
| 
Im 
nal 


ita 


E 
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3.2 — Matriz de deformabilidade 
3.2.1 — Matriz de deformabilidade do terreno 


Consideremos a carga triangular de resultante 


unitária representada na fig. 3 c). O elemento E 
da matriz de deformabilidade do terreno é o deslo- 
camento produzido por esta carga unitária distribuida 
entre os nós (j-1) e (j4y 1), no nó 1. A densidade 
P, será portanto: 


1 


= —— (3.2.1.1) 
a(Z;p— 2-1) 


P; 


O deslocamento fi será então calculado por so- 
breposição dos efeitos correspondentes às cargas (1) 
e (O indicadas, a partir da equação (2.2.3). As car- 
gas (T) e (2) ficam definidas por: 


pO=—(p Z;-)Mz;- spo (P; Zi) [(Zn— 2) 


D D 
PTB KZ 21) po =-p [(Z1— 2) 
(3.2.1.2) 
pelo que: 


t D D 
f=u(Z v2 DP! Z, ) + 


7) 
pe u(z; ' Zu , po , ZA , Hi ) (3.2.1.3) 


A matriz de deformabilidade do terreno será de- 
signada por: 


P=ltl o 4)=14,8 un (3.2.1.4) 


3.2.2 — Matriz de deformabilidade da estaca 


Consideremos a carga triangular unitária repre- 
sentada na fig. 3 d). O deslocamento do nó i pode de- 
compor-se em duas parcelas. A primeira corresponde 
a um deslocamento de corpo rígido da estaca devido 


à deformabilidade dos «apoios» extremos É pp À se 


f 
gunda corresponde à flexão da estaca f, = 
Quanto à primeira parcela temos: 


f — f 
) depan dy | (3.2.2.1) 


em que A ef são os deslocamentos no terreno 


| (n+1)3 
nos nós o e (n+1) produzidos por cargas iguais mas 
de sentidos contrários às reacções distribuidas, 
fig. 3 d). 


Deste modo f, MM: 


; (n+0) serão calculadas por 
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(2.2.2). As densidades Es e F(n+1)) determinadas de 
modo a equilibrar a carga triangular unitária são: 


I 1 
Toj = Po éra — Foi) ( Z,1— Zy | (3.2.2.2) 


E (o 41) = [— ZP(Z1— 2% 1)+P; [(2 — 2% MZ)-1 + 
+ + (2 — 23 22 +) Lapa Za MZa + 
+22) —U(Z041— 2)) (3.2.2.3) 


Teremos então de acordo com (2.2.2) e (2.2.2.3): 


r oj 
É; a u(0,Z,,Foj E E , 0) + u (Zu ] l , 
| 


Ta + 1)jZa 


r 
(041) 5 
E 


(3.2.2.4) 
|— z, 


ro) 
Eat) = Uu (9,255; 8: E 1) + u (Za , , 


Tn + 1)j Zn 
l- z, 


F(n+41)j 


ca 1) (3.2.2.5) 


A segunda parcela a corresponde à flexão da 
estaca e representa portanto o deslocamento produ- 
zido no nó i por uma carga triangular unitária e 
pelas correspondentes «reacções» distribuidas de 
apoio. As condições de fronteira correspondem à viga 
simplesmente apoiada em z=o e z=1 ou seja 


d'u du 


O elemento £f,, da matriz de deformabilidade da es- 
taca será então: 


e a f 
to, =fy +%, (3.2.2.6) 
A matriz é designada por: 
Pe=[f,] 14)=128,..,n (3227) 


3.2.3 — Compatibilidade 


A compatibilidade dos deslocamentos do solo e 
da estaca no nó j exprime-se por: 


n 
t a 
j=1 
(3.2.2.8) 
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em que a componente Q, do vector das incógnitas, 
Q = (Q;), representa o valor da resultante da carga 


triangular distribuída entre os nós (j-1) e (j+1). 
A equação (3.2.2.8) pode tomar a forma matri- 

cial: 
(Pr +Fe) Q=-=U, (3.2.2.9) 


3.2.4 — Cálculo dos deslocamentos e esforços na 
estaca 


Atendendo à compatibilidade, os deslocamentos 
finais na estaca podem ser facilmente obtidos uma 
vez conhecido o vector Q, pela equação (3.2.2.9). 


ue = ut = FtQ (3.2.4.1) 


em que ue e ut são os vectores dos deslocamentos 


nodais na estaca e no solo, respectivamente. O cál- 
culo dos esforços na estaca pode ser feito directa- 
mente a partir do vector Q, ou a partir dos desloca- 


mentos u*, sendo em qualquer caso um simples pro- 


blema de Resistência de Materiais, pelo que aqui não 
será abordado. 


4 — PROGRAMA DE CALCULO AUTOMÁTICO — APME 
— Aplicação do Problema de Mindlin a Estacas 


Compõe-se fundamentalmente dum programa 
principal e quatro subrotinas. Duas delas, «SUBROU- 
TINE UTE> e <«<SUBROUTINE DELTA», calculam 
as matrizes de deformabilidade do terreno, F'*, da es- 


taca, F*, e o vector dos deslocamentos no sistema base, 
U., tendo em conta as equações expostas em (2.2), 


(3.1) e (3.2). As duas restantes, <SUBROUTINE 
ARRAY» e «SUBROUTINE SIMQ» do CCUTL(*), 
destinam-se simplesmente à resolução do sistema de 
equações lineares (3.2.2.9). 

O programa APME foi desenvolvido em lingua- 
gem FORTRAN IV para o Computador IBM 360/DOS 
do CCUTL, O número total de memórias ocupadas 
é de 51 K bytes, para um número máximo de 30 nós 
intermédios. 

O tempo de execução do programa para os exem- 
plos referidos nunca ultrapassou os 3 minutos, Apre- 
senta-se em anexo um organograma do programa 
APME. 


5 — EXEMPLOS 


Apresentam-se em anexo os diagramas de des- 
locamentos em estacas obtidos com o programa 
APME, os quais são comparados com os obtidos em 
[2] pelo método dos elementos finitos. Refere-se que 
estes últimos foram obtidos para estacas de secção 
circular pelo que, tendo sido consideradas neste traba- 


(*) Centro de Cálculo da Universidade Técnica de Lisboa. 
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lho estacas rectangulares, se tomou para efeitos de 
comparação uma estaca rectangular de largura igual 
ao diâmetro e com o mesmo momento de inércia. 


6 — CONCLUSÕES 


Pretendemos com esta comunicação apresentar 
um método capaz de permitir a comprovação dos re- 
sultados obtidos pela técnica dos elementos finitos 
na resolução do problema em causa, nos casos em 
que o meio é homogéneo, resistente à tracção e de 
comportamento linear, Esta comprovação pareceu-nos 
necessária por julgarmos que o número de elementos 
empregados é muitas vezes insuficiente. E, efectiva- 
mente, os resultados obtidos mostram certas discre- 
pâncias entre os resultados obtidos pelos dois méto- 
dos. 

O método que apresentámos pode servir de afe- 
ridor porque o erro de discretização com ele asso- 
ciado é geralmente muito inferior ao correspondente 
ao método dos elementos finitos. Lembramos que o 
referido erro provém, no nosso caso, do facto de a 
compatibilização dos deslocamentos ser feita exclu- 
sivamente na direcção da força aplicada na cabeça 
da estaca e somente num número finito de nós ao 
longo do respectivo eixo. Esta última causa de erro 
pode porém ser facilmente controlada fazendo variar 
o número de nós e observando as consequências desta 
variação. 

Está evidentemente fora de questão a grande 
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generalidade do método dos elementos finitos, que lhe 
dá possibilidades de simulação do comportamento real 
da estaca e do terreno com que o nosso não pode 
competir. Não queremos no entanto deixar de chamar 
a atenção para a não resistência à tracção do solo 
poder ser facilmente simulada com a nossa técnica 
utilizando o método simplificado proposto por Ka- 
nayan (ver [2]). 
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Com base num modelo apropriado dum dispositivo reten- 
tor de amostras, estabelece-se a expressão da função densi- 
dade espectral à saida, devida à presença à entrada, de flu- 
tuações com função de auto-correlação especificada. 

O modelo utilizado para representar o retentor leva em 
linha de conta os valores não desprezáveis da constante de 
tempo e do tempo de captação do dispositivo. 

Apresentam-se os resultados de medidas efectuadas num 
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deduzem por integração os deslocamentos devidos a uma 
carga distribuida linearmente, apresenta-se um método para 
a resolução do problema da estaca cravada num solo homo- 
géneo solicitada por uma força horizontal. 
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gona! coordinates are used and the method can be mecha- 
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on the different coordinate planes. 
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title «Geometrical Method to obtain Lorentz Transformations. 
Properties of the Method and its Applications». 
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Based on an appropriate model of a sample-and-hold 
device, a theory is developed leading to an expression for 
the output spectral density function due to the presence at 
the input of fluctuations with a specified auto-correlation 
function 

The sample-and-hold mode! used takes into account the 
non-negligible time constant and acquisition time of the 
device, 

Results of measurements made on a real device are 
presented to confirm the theoretical predictions and an appli- 
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